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ABSTRACT 


The  use  of  aquatic  macroinvertebrates  for  monitoring  water  quality  has  become 
popular.  This  paper  assessed  the  ability  of  rapid  assessment  using  a one  minute  riffle 
kick  net  method  to  detect  natural  within-stream  and  seasonal  macroinvertebrate 
variability  in  the  mountain  stream  ecoregion  of  Montana.  It  also  assessed  the  ability  of 
these  methods  to  detect  water  quality  impairment.  Results  suggest  that  the  methods 
used  are  able  to  reflect  natural  variability  and  that  some  indices  or  'metrics'  are 
potentially  more  useful  than  others  for  detecting  water  quality  impairment  in  heavily 
•impaired'  streams.  It  is  unkown  how  useful  these  methods  are  for  detecting  water 
quality  impairment  in  less  severely  impaired  streams.  Suggestions  are  also  made  on 
how  to  increase  the  sensitivity  of  these  methods. 
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INTRODUCTION 

Freshwater  is  one  of  the  most  important  resources  in  the  world,  however, 
within  the  last  two  hundred  years,  technology  and  human  population  growth  have 
caused  tremendous  impacts  on  freshwater  ecosystems  (Bau  1995;  Karr  1991).  For 
example,  Benke  (1990)  reported  that  98%  of  the  3.2  million  miles  of  rivers  in  the  U.  S. 
were  not  healthy  enough  to  be  considered  high  quality  and  worthy  of  federal 
protection. 

The  U.S.  Clean  Water  Act  of  1972  (CWA)  reflects  the  concerns  of  the  people 
of  the  United  States  regarding  clean  water  and  healthy  aquatic  ecosystems  and 
mandates,  as  its  goal,  "...to  improve  and  maintain  the  physical,  chemical,  and 
biological  integrity  of  our  nation's  waters".  Biological  integrity  can  be  defined  as 
"...the  ability  of  an  aquatic  ecosystem,  to  support  and  maintain  a balanced,  integrated, 
adaptive  community  of  organisms  having  a species  composition,  diversity,  and 
functional  organization  comparable  to  that  of  the  natural  habitats  of  a region"  (Karr  and 
Dudley  1981). 

The  U.  S.  Environmental  Protection  Agency  (U.S.  EPA)  is  the  federal 
regulatory  agency  responsible  for  the  enforcement  of  the  Clean  Water  Act  of  1972.  The 
U.S.  EPA  has  encouraged  individual  states  to  develop  their  own  regionalized  methods 
of  assessing  the  biological  integrity  of  their  (states)  surface  waters  based  on  a model 
proposed  in  the  U.S.  EPA's  Rapid  Bioassessment  Protocols  (RBP's)(U.S.  EPA  1989). 


? iv:  7 


^ . ' i 


■ '■  V .m 


,A> 


l^on'^uQp^T..; 


;■) 

■\u4  iktifi'tfx  W!*;  V;  ^ • •■  '/  ;V».  .<ft^  r^^iiiy/;;  /V 

“ ..  .'  ■ .. 


“ * * • 

' ' * '''  , ‘ 

’>r?<  fi'-.  I ) 

• , -.v:  .1 

- ' "V; 

nfi,  i<-  t‘jv.ix  ;u  •:  :4]  . -'our.  V ■■»■  ‘ij  j"..  B.U 

' ’ " > ‘'f  ’ flaSUfj  ' ^ ■ 


T i ''’'  ^ iCKT'  'wj'  r5,'- 

i -••  ■•  ■':'  I'v  ..  'ir'  ' Ai  '■  ■'  ’’  ' 

u'»  iv) '.  ,\JK>w  ,r  r '‘  i-  it%  *.■'■.  •>'.:  ■ 


u<:‘  'is.ViT'.’tfi.*  »v  i.4j|)r  V{ji'^“..'  ',  .'j  *!''■■;'•  .•  :^if:W.r,^Ai  lo- 


{ins.  I'M'.'  «.'  I/,-  /.'-v**?**  ••'l  ti  - /I  K^\.x  o.'  . * •,'•#«•:{  >ji  ?&  . 

\ ^ W#«  ' 


’■-i,\  ’ ■,.  ’ 'I  i.i'.r;  ih  Mi: 


9 ” i-’j.  *>  .'  . rvri-j  ' {?/■■■•'•'  vlmvirn 


•i. 

^ '-■• ' V'^v  ■ :; 
"U-  ' 

.if*  ?f  \i\‘yd  .?  *")  vp?Dj|A  t 

I 1 '• 

:<J‘r  i ^ jv  -s.'-Lv  lo'i 

' ■ !■  '•  J " ' 


iiwni'i^n  faii'i-'uoi;.,m  ir  /n  u;r'  ikyUnt)  i-iWp  '.iulihitini  U’{ 


j»>:r Jg;,:T 


, ) . , .,|... 

, >'^'W  ■ >Hi  . '.v*-'’-'  -Jiu*  iMU  V)  feulgyk  ’^*  'jijj 

t 

-^s»i»n’rA^  ./’'I' , i Mi  iTH.  i i Ui:^ti.^p.r;iL)50^id  *h:i 


2 


Most  RBP’s  rely  on  a multi-metric  approach  using  reference  conditions  from  a 
regional  framework  (Davis  and  Simon  1995).  A metric  is  simply  a numeric  index 
which  estimates  some  attribute  of  a community,  for  example,  the  number  of  taxa 
present  in  a sample  or  the  number  of  individual  organisms  in  a sample.  Metrics  that 
were  used  in  this  study  are  in  Table  1. 


Table  1.  Metrics  used  in  this  study. 


Taxa  richness* 

EPT  Richness'* 

% dominant  taxa 

it  of  organisms/  sample 

% EPT  individuals 

modified  HBP 

it  of  Ephemeroptera  taxa 

QSI  for  taxa^ 

Shannon  Index^ 

it  of  Plecoptera  taxa 

it  of  eg  and  fg  taxa^ 

% predator  taxa 

it  of  Trichoptera  taxa 

it  of  shredder  taxa 

% shredder  taxa 

% semivoltine  taxa 

a of  scraper  taxa 

% scraper  taxa 

% multivoltine  taxa 
% univoltine  taxa 

it  of  predator  taxa 
MTr 

% eg  and  fg  taxa 

‘Taxa  richness  is  the  number  of  different  taxa 

^ EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa 

^ QSI  for  taxa  is  the  Quantitative  Similarity  Index  based  on  percent  relative  abundances  of  taxa  between 

samples 

^ # of  eg  and  fg  taxa  is  the  number  of  collector-gatherer  and  filterer-gatherer  taxa  based  on  Merrit  and 
Cummins  (1984) 

^ HBI  is  the  Hilsenhoff  Biotic  Index  based  on  organic  pollution  tolerance  values  (Hilsenhoff  1987  and 
Bukantis  1995b) 

^ Shannon  Index  is  a biotic  diversity  index 
‘‘MTI  is  a metals  tolerance  index  (Bukantis  1995b) 


Reference  biological  and  habitat  conditions  are  obtained  from  "least  impaired" 
sites  in  an  ecoregion  or  sub-ecoregion  (regional  framework)  as  defined  by  Omemik 
(1987).  A biological  assessment  is  then  made  for  the  site  that  is  thought  to  be  impaired 
within  the  reference  ecoregion  using  this  multi-metric  approach.  RBP’s  are  used  by 
many  state  and  federal  agencies  for  monitoring  water  quality.  Regionalized  RBP's  are 
designed  to  be  efficient,  economical,  easily  interpretable,  and  accurate  in  assessing 
water  quality  impairment,  at  least  on  a gross  level  (Resh  and  Jackson  1993).  Results  of 
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many  RBP's  can  be  implemented  in  less  than  5 days  (Lenat  and  Barbour  1994) 
although,  often  this  may  not  be  true  for  agencies  that  do  not  have  the  resources 
available  (Bukantis  1995).  RBP's  can  be  used  for  both  point  and  non-point  sources  of 
impairment  (Resh  and  Jackson  1993).  Aquatic  macroinvertebrates  are  the  most 
frequently  used  organism  in  RBP  assessments  (Rosenberg  and  Resh  1993). 

The  State  of  Montana  is  presently  developing  its  own  modified  RBP's  for 
wadeable  rivers  and  streams  in  Montana  and  has  delineated  the  state  into  three 
ecoregions:  1)  the  plains,  2)  foothills  and  intermountain  valleys,  and  3)  mountain 
ecoregions  based  on  Omemik  (1987)  and  from  preliminary  findings  of 
macroinvertebrate  communities  in  Montana  using  Detrended  Correspondence  Analysis 
(DECORANA)  (Wisseman  1990).  Water  quality  impairment  in  these  three  ecoregions 
in  Montana  is  mostly  due  to  non-point  sources  (Montana  305b  Report  1994).  The 
majority  of  non-point  impairment  in  Montana  to  rivers  and  streams  is  from  agriculture, 
logging  and  mining  (particularly  in  the  mountain  ecoregion),  road  building, 
sedimentation,  irrigation,  and  habitat  loss  (Montana  305b  Report  1994).  Of  the 
176,750  miles  of  streams  in  Montana,  only  10%  have  been  assessed  for  water  quality 
(Montana  305b  Report  1994).  The  State  of  Montana  Water  Quality  Division  of  the 
Department  of  Environmental  Quality  has  been  developing  RBP's  since  1989 
(Wisseman  1990)  and  is  increasing  its  data  base  annually  (McGuirre  1994,  1995). 
Bollman  (1995)  is  presently  completing  an  extensive  study  in  the  foothill  and 
intermountain  ecoregion  in  Montana.  Very  little  data  exists  using  macroinvertebrates 
to  assess  water  quality  impairment  in  the  mountain  streams  ecoregion  of  Montana. 
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HYPOTHESES  TESTED 
I tested  four  hypotheses; 

1)  within-site  variation  exists  in  the  abundance  and  distribution  of 
macroinvertebrate  taxa  in  mountain  streams  of  Montana  and  can  be 
detected  with  methods  used  in  this  study. 

2)  within-stream  variation  exists  in  the  abundance  and  distribution  of 
macroinvertebrate  taxa  in  mountain  streams  of  Montana  and  can  be 
detected  with  methods  used  in  this  study. 

3)  seasonal  variation  exists  in  the  abundance  and  distribution  of 
macroinvertebrate  taxa  in  mountain  streams  of  Montana  and  can  be 
detected  with  methods  used  in  this  study. 

4)  water  quality  impairment  affects  the  abundance  and  distribution  of 
macroinvertebrate  taxa  in  mountain  streams  of  Montana  and  can  be 


detected  with  methods  used  in  this  study. 
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LITERATURE  REVIEW 
Water  Quality  Problems  in  the  U.S. 

Human  caused  degradation  of  river  ecosystems  worldwide  has  reached 
unprecedented  levels.  "Not  one  riverine  system  in  America  has  been  spared "(Doppelt 
et  al.  1993).  Dopplet  et  al.  (1993)  report  that  50%  of  our  nations'  waters  fail  to  meet 
federal  water  quality  standards.  Also,  less  than  2%  of  our  nations 's  rivers  even  qualify 
for  Wild  and  Scenic  designation.  From  one-third  to  three-fourths  of  aquatic  species 
nationwide  are  rare  to  extinct.  Aquatic  species  are  disappearing  at  a faster  rate  than 
terrestrial  species  (Doppelt  et  al.  1993)  with  twenty  percent  of  native  fishes  in  the 
western  U.S.  extinct  or  endangered  (Miller  et  al.  1989).  Two  hundred  and  fourteen 
stocks  of  native  Pacific  salmon  and  steelhead  are  facing  extinction  (Nehlsen  et  al.  1991) 
and  commercial  fish  harvests  in  major  rivers  throughout  the  U.S.  have  declined  from 
80%  to  100%  (Karr  et  al.  1985;  Ebel  et  al.  1989;  and  Patrick  1992).  Karr  (1991) 
asked  'what  would  our  country  think  if  U.  S.  agricultural  productivity  decreased  by 
80%?' 

It  is  estimated  that  60%  to  80%  of  natural  riparian  vegetation  has  already  been 
lost  or  is  degraded  due  to  human  activities  in  the  U.S.  (Swift  1984),  while  seventy 
percent  of  our  nations'  rivers  and  streams  have  been  impaired  by  flow  alteration 
(Doppelt  et  al.  1993).  Over  600, CKX)  miles  of  rivers  have  already  been  dammed  in  the 
U.  S.,  with  the  Yellowstone  River  being  the  only  remaining  free-  flowing  large  river 
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(longer  than  6(X)  miles)  left  in  the  U.S.,  outside  of  Alaska  (Karr  1995). 

The  only  attempt  at  quantifying  river  channelization  in  the  U.S.  was  conducted 
over  20  years  ago  in  1972  by  the  Council  of  Environmental  Quality  which  estimated 
that  235,000  miles  of  streams  and  rivers  in  the  U.S.  have  been  channelized  (Doppelt  et 
al.  1993).  The  United  States  Environmental  Protection  Agency  (U.S.  EPA)  reports  that 
nonpoint  sources  affect  65%  of  impaired  streams  (Karr  1991). 

Water  Quality  Problems  in  Montana 

Assessments  of  water  quality  in  Montana  have  been  completed  in  only  10%  of 
its  stream  miles  ( Montana  Water  Quality  Division  1994).  The  State  of  Montana 
classifies  its  streams  as  fully  supporting,  partially  supporting,  or  non-supporting  of 
their  designated  uses.  Non-supporting  stream  miles  total  925  (about  5%  of  those 
assessed)  and  partially  supporting  streams  total  75  % of  those  assessed.  It  is  unknown 
how  many  unassessed  stream  miles  (159,070  miles)  are  fully  supporting  their 
designated  uses.  Agriculture  has  impaired  75%  of  the  stream  miles  that  have  been 
assessed  and  90%  of  stream  mileage  impairment  in  Montana  comes  from  non-point 
sources  (Montana  Water  Quality  Division  1994).  Nutrients,  siltation,  suspended 
solids,  salinity,  flow  and  habitat  alterations,  and  metals  are  the  major  causes  of  stream 
impairment  in  Montana  ( Montana  Water  Quality  Division  1994). 
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Biomonitoring  Approaches 

The  idea  of  biological  indicator  organisms  is  not  new.  A kings'  use  of 
winetasters  and  the  use  of  canaries  by  miners  to  detect  dangerous  gases  are  familiar 
examples  of  biological  indicators  (Cairns  and  Pratt  1993).  Even  the  relative  abundance 
of  rats  in  our  cities  is  a biological  indication  of  environmental  conditions. 

The  use  of  freshwater  aquatic  organisms  as  bioindicators  of  water  quality 
appears  to  have  its  origins  in  Europe  with  the  work  of  Kolenati  in  1848.  Kolenati 
(1848)  reported  that  city  effluents  caused  the  disappearance  of  caddisflies  downstream. 
Kolkwitz  and  Marsson  (1908)  related  the  degree  of  pollution  in  a river  as  a measure  of 
sewage  organic  matter.  Other  researchers  then  began  developing  lists  of  pollution 
'indicator  organisms'  in  rivers  and  streams  (Richardson  1925,  1929;  Gaufm  1958; 
Cairns  and  Pratt  1993).  Cairns  and  Pratt  (1993)  suggest  that  presently,  "biological 
surveillance  of  communities  - with  special  emphasis  on  characterizing  taxonomic 
richness  and  composition  - is  perhaps  the  most  sensitive  tool  now  available  for  quickly 
and  accurately  detecting  alterations  in  aquatic  ecosystems". 

Aquatic  macroinvertebrates  are  the  most  commonly  used  indicator  organisms 
for  water  quality  assessment  (Rosenberg  and  Resh  1993),  although  other  aquatic 
organisms  are  routinely  used,  including  fishes  (Simon  and  Lyons  1995)  and  periphyton 
(Bahls  1993).  Some  of  the  reasons  for  using  aquatic  macroinvertebrates  as  indicators  of 
water  quality  include;  1)  macroinvertebrates  occur  in  all  types  of  freshwater 
environments  (Merrit  and  Cummins  1984),  2)  the  many  different  kinds  of 


a ,^3n  5m  'Kji4:i,‘^#..t^ 


uiiiUnu.'’  '''.■  <!'  ::x^n 

•i;.«  •►'  ^iViU  i'^btUn  r,iiin-K^  ,{t^{,'  f:4’i«i  h^,n  7:m'iA 

1 - j , ■ 


i.sk&s-lv 


, , _ ,,, 

m-L-Uthmo  i*s:r-*='U)W»/«v/<3.‘io  tfu-rmnm' 

I ,.  : ^ 'ji^W^*^mm 


■ ■ . ■ . .'  . V''t'-^^^  i 

■.  : / ' ■ ' 

■ ' ‘ • ' ' ' ' ' ’ir''  ' 


■I 


'fe' 


^ xai:, 

jTO  ' 


Vr' 


'X.''* ■>/ ici <iil  , v^)[ ‘ rjfs<u'?^.^  jfie'^rr7, a 


J. 

•>:■<■  4*-. 


n r.iffii  vj  . 'Woe  imf  -w JiH  ..r^#C3^4<^-2!^ 

. ■ • ' ^ •'•  i . T-- 


' 


&>ii\ 


IT-:  Ite  X{i|r>*ii*b  $■: 

^ ^ '^.|:3':&  T-54, 


■mi  'hi  iMBTjwt  Hj  ni  yitxia  atnf  '^4‘tiwMrmn  < i ■^|fcoi^fe 


\ 1'; 


m ’ \V 

jjQ  u.  i ■\^ . 


M. 


■Jt>  ItiJOh  M*'iA>>.vi)  . 'lev'll  (^C 

m.  I ' . ' ■ w 


■ " w •/ 

■i  .,§.l(W’-i 


V.',  y) 


A,  .1 


A '■< } ^ 

Ji  t 


i5y 

, :t;^5i  . 4 

Ikiii^ 


8 


macroinvertebrate  species  found  in  freshwater  habitats  allows  for  a wide  range  of 
responses  to  different  sources  of  impairment  (Hellawell  1986,  Abel  1989),  3)  aquatic 
macroinvertebrates  are  less  mobile  than  fishes  which  allows  for  spatial  analysis  of 
impairment  (Slack  et  al.  1973,  Hellawell  1986)  and  4)  macroinvertebrates  have  longer 
life  cycles  compared  with  microorganisms  and  periphyton,  which  allows  for  analysis  of 
disturbance  over  time  (Gaufm  1973,  Abel  1989).  Aquatic  macroinvertebrates  are  "on 
site  monitoring  the  water,  regardless  of  the  presence  or  absence  of  the  investigating 
biologist"  (Nehrig  1976). 

Many  types  of  biomonitoring  of  freshwater  ecosystems  are  in  use  (Johnson  et. 
al  1993).  The  use  of  individual  organisms  for  toxicity  testing  for  various  pollutants  in 
the  laboratory,  in  artificial  streams,  and  in  field  studies  is  common  (Gaufm  1973,  Rahel 
and  Kolar  1990,  Dunkel  and  Richards,  submitted).  These  studies  typically  measure 
physiological,  biochemical,  behavioral,  and  life  history  changes  of  the  individual  test 
organisms  to  various  types  of  pollutants  (Marten  and  Zwick  1989,  Nardi  and  Watson 
1990,  Johnson  et  al.  1993). 

Population  and  community  structure  studies  are  also  used  to  measure  the  effects 
of  different  pollutants  on  aquatic  macroinvertebrates  both  in  the  laboratory  and  more 
commonly  in  field  studies.  These  studies  may  use  a univariate  or  multivariate  approach 
or  a combination  of  the  two  (Johnson  et  al.  1993). 

Univariate  studies  of  populations  and  communities  can  be  used  to  develop  biotic 
indices  and  scoring  systems.  An  example  of  such  an  index  is  the  Hilsenhoff  Biotic 
Index  which  scores  the  sensitivities  of  many  aquatic  macroinvertebrates  to  organic 
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pollution  and  is  widely  used  in  monitoring  programs  (Hilsenhoff  1987,  Resh  and 
Jackson  1993).  A prerequisite  to  the  use  of  biotic  indices  is  a thorough  knowledge  of 
how  the  type  of  pollution  used  for  each  index  affects  the  biota  (Johnson  et  al  1993). 

Often  graphical  analysis  of  macroinvertebrate  data  can  provide  more  insight  into 
biology  than  a simple  p-value  (Fore  et  al.  1995).  For  example.  Fore  et  al.  (1995)  were 
interested  in  whether  scores  (metrics)  for  best  and  worst  sites  clumped  together  tightly 
or  were  spread  more  evenly.  Graphs  allowed  them  to  examine  each  metric's  range  and 
evaluate  where  along  the  continuum  of  degradation  the  metric  was  most  sensitive.  Fore 
et  al.  (1995)  suggest  that  the  main  problem  with  statistical  correlation  is  that  no  single 
variable  can  summarize  all  the  human  activities  that  degrade  streams.  Different  data 
may  be  appropriate  for  different  types  of  degradation.  Fore  et  al.  (1995)  doubt  that 
any  single  land-use  or  chemical  variable  can  reliably  rank  sites.  Fore  et  al.  (1995) 
suggest  that  this  is  a primary  reason  for  evaluating  the  river  condition  according  to  the 
resident  biota.  In  most  cases,  sites  that  appear  to  a human  observer  to  be  heavily  used 
and  visibly  damaged  usually  are  more  degraded  than  those  that  appear  pristine  (Fore  et 
al.  1995). 

Multivariate  approaches  including  regression  analysis  can  be  used  to  ordinate 
populations  or  communities  along  environmental  gradients  (Digby  and  Kempton  1994). 
Johnson  et  al.  (1993)  suggested  that  direct  gradient  analysis  (regression)  could  be  used 
to  predict  species  assemblages  along  a known  environmental  gradient.  This 
information  can  then  later  be  used  to  infer  water  quality  from  species  assemblages 


found  in  a study  site. 
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Principal  Components  Analysis  (PCA)  is  an  example  of  indirect  gradient 
analyses  which  measure  species  assemblages  gradients  regardless  of  any  environmental 
gradient  (Johnsen  et  al.  1993).  For  example,  in  Principal  Components  Analysis, 
Principal  Component  1 (or  eigenvector  1)  contains  the  largest  individual  differences  in 
a species  assemblage.  The  second  Principal  Component  (eigenvector)  accounts  for  the 
largest  proportion  of  remaining  individual  differences  not  correlated  to  (or  orthogonal 
to)  Principal  Component  1 (Williams  and  Feltmate  1992).  PCA  calculates  the  line,  or 
component,  that  extracts  the  maximum  amount  of  statistical  variance  from  a cloud  of 
points  (Tabachnick  and  Fidell  1989).  For  stream  macroinvertebrate  analysis,  each 
point  then  will  represent  a stream  site.  The  number  of  dimensions  through  which  the 
line  passes  would  equal  the  number  of  taxa  collected  (Fore  et  al.l995).  Usually,  PCA 
data  uses  species  lists  and  abundances  to  interpret  differences  between  stream  sites 
(Norris  and  Georges  1993).  If  land  use  data  are  not  sufficient  to  perform  canonical 
correspondence  analysis,  then  PCA  may  be  used  and  site  characteristics  can  later  be 
identified  on  the  principal  components  plots  (Fore  et  al.  1995).  Hannaford  and  Resh 
(1995)  used  PCA  to  analyze  the  correlation  matrix  of  six  RBP  metrics  used  in  their 
study  of  variability  in  rapid-bioassessment  surveys.  In  an  intensive  ecological  study, 
Gustafson  (1990)  was  able  to  use  principal  components  analysis  to  identify  the  velocity 
preferences  of  aquatic  insects  in  the  Gallatin  River,  MT  and  identify  the  insect  species 
distribution  to  three  distinct  sections  of  the  river.  Gustafson  (1990)  also  used  a faunal 
and  physical  component  PCA  analysis  to  easily  identify  headwater  insect  species,  lower 
Gallatin  River  species,  and  widespread  species  within  the  drainage. 
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With  the  many  types  of  data  analysis  available,  it  should  be  pointed  out  that 
..."Too  often  biological  patterns  in  complex  data  sets  are  reported  in  terms  of 
complicated,  high  level  statistics  such  as  principal  components  analysis"  (Fore  et  al. 
1995).  While  many  times  it  is  easier  to  compare  sites  and  search  for  patterns  in  the  raw 
data  by  simply  noticing  which  taxa  were  present  or  absent  in  similar  sites  (Fore  et  al. 
1995).  All  of  these  quantitative  methods  mentioned  above  require  many  replicates  with 
detailed  statistical  analysis,  which  may  be  cost  prohibitive  given  the  thousands  of  miles 
of  rivers  and  streams  potentially  impaired  in  the  U.  S.  (Resh  and  Jackson  1993). 

Rapid  Assessment 

With  the  high  cost  associated  with  quantitative  methods,  limited  budgets  and 
limited  personnel  in  management  agencies  responsible  for  water  quality,  and  given  the 
increasing  deterioration  of  water  quality  throughout  the  world,  the  most  recent 
approach  in  biomonitoring  has  been  "rapid  assessment" . Rapid  assessment  is  a 
qualitative  method  (sometimes  semi-quantitative)  designed  to  measure  water  quality  and 
can  also  be  used  to  measure  long-term  regional  changes  (Resh  and  Jackson  1993). 

63%  of  the  123  U.S.D.A.  National  Forest  Districts  and  61%  of  the  56 
Department  of  Interior,  Bureau  of  Land  Management  Districts  surveyed  by  Angradi 
and  Vinson  1995,  reported  conducting  some  type  of  aquatic  macroinvertebrate 
monitoring.  Mining,  timber  harvest,  and  grazing  management  were  the  most  often 
cited  reasons  for  using  aquatic  macroinvertebrate  monitoring.  Protocol  IE  is  the  most 
rigorous  tier  level  approach  for  the  RBP's  suggested  by  the  U.S.  EPA  and  is  used  by 
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over  20  states  (Resh  et  al.  1995).  Protocol  HI  identifies  macroinvertebrates  to  the  genus 
or  species  level  and  incorporates  a habitat  assessment  for  each  site. 

Many  of  the  assessment  methods  used  throughout  the  U.S.  rely  on  a multimetric 
approach,  while  in  Europe  a multivariate  approach  is  more  common.  Norris  (1995) 
discussed  the  use  of  multivariate  methods  in  Europe  and  suggested  that  multivariate 
methods  are  seriously  under  used  in  the  U.S.  Gerritsen  (1995)  on  the  other  hand, 
argues  that  multimetric  methods  used  in  the  U.S.  function  well  and  that  multivariate 
methods  "..are  more  complex,  require  additional  specialization  by  practitioners,  and 
are  more  difficult  to  convey  to  managers  and  the  public".  Fore  et  al.  (1995)  state  that, 
"multimetric  indices  formalize  what  any  good  biologist  familiar  with  local  biota,  knows 
about  the  biological  condition  of  a stream".  Indices,  thus,  become  important  tools  for 
communication  with  non  specialists.  A multimetric  approach  can  be  compared  to 
familiar  economic  indexes  such  as  the  consumer  price  index  which  has  been  widely 
used  for  over  50  years  (Fore  et  al.  1995). 

Rapid  bioassessments  are  frequently  compared  to  the  use  of  thermometers  in 
assessing  human  health  (Karr  1995;  Resh  and  Jackson  1993).  Results  are  easily 
obtained  and  can  be  compared  to  the  "normal"  condition.  Resh  and  Jackson  (1993) 
point  out  that  the  key  questions  then  become:  what  populations  and  community 
measures  are  biologically  relevant  (the  thermometers)?,  what  are  the  thresholds  against 
which  they  are  being  compared  (the  normal  body  temperature)?  and,  how  much  of  a 
deviation  from  the  threshold  is  a sign  of  "ill  health"?  They  suggest  that  a critical  and 
often  controversial  decision  is  the  selection  of  the  most  appropriate  measures  for  use  in 
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multimetric  rapid  assessment  approaches. 

Many  metrics  or  indices  are  presently  being  tested  and  used  in  rapid  assessment 
programs  throughout  the  U.S.  including  the  Pacific  Northwest  (Fore  et  al.  1995).  Resh 
and  Jackson  (1993)  gave  an  excellent  overview  of  the  many  metrics  that  have  been 
proposed  for  use,  as  well  as  the  various  collection  methods,  sampling  strategies  and 
taxonomic  resolution  used  by  various  researchers  and  agencies. 

Problems  with  Rapid  Assessments 

Even  though  rapid  assessments  are  becoming  very  popular,  their  ability  to 
detect  water  quality  impairment  is  often  questioned  (Norris  1995,  Southerland  and 
Stribling  1995,  Resh  and  Jackson  1993).  Although  rapid  bioassessments  were  designed 
to  be  cost  effective,  several  Forest  Service  Districts  and  Bureau  of  Land  Management 
Districts  surveyed  by  Angradi  and  Vinson  (1995)  reported  that  they  could  not  afford  to 
conduct  assessments  properly.  Several  districts  also  reported  that  natural  variability  is 
so  great  in  aquatic  macroinvertebrate  communities  that  current  methods  are  unreliable 
or  insensitive  to  land  management  impacts. 

Another  problem  using  RBP's  concerning  natural  within-site  variability  is  that 
most  RBP's,  including  Protocol  El  of  the  U.S.  EPA,  suggest  collection  of  only  one 
sample  per  site  (Resh  and  Rosenberg  1989).  Research  has  shown  variability  in  the 
replication  of  many  RBP  metrics  (Resh  et  al.  1990,  Steven  and  Szczytko  1990,  Barbour 
et  al.  1992,  Hannaford  and  Resh  1995,  and  Resh  1994). 

The  choice  of  sampling  methods  and  the  choice  of  the  most  appropriate  habitat 
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to  be  sampled  may  influence  results  (Rosenberg  and  Resh  1993).  Even  the  choice  of 
the  portion  of  a riffle  to  be  sampled  may  be  biased  due  to  the  distribution  of  benthic 
organisms  within  the  riffle.  Brown  and  Brown  (1984)  showed  a strong  upstream  biased 
distribution  of  insects  within  riffles.  With  all  of  the  problems  seemingly  associated  with 
rapid  assessment,  many  researchers  and  mangers  feel  that  these  problems  can  be 
addressed  and  evenmally  overcome  with  an  increased  understanding  of 
macroinvertebrate  community  responses  to  water  quality  impairment. 

Rapid  Assessment  in  Montana 

The  Water  Quality  Division  (WQD)  of  the  Montana  Department  of 
Environmental  Quality  (DEQ)  is  presently  assessing  the  applicability  of  its  modified 
version  of  the  U.  S.  Environmental  Protection  Agencies  (EPA)  Rapid  Bioassessment 
Protocols  (RBP’s)  using  macroinvertebrates  for  monitoring  water  quality  in  mountain 
streams  of  Montana.  These  methods  use  a multi-metric  approach  based  on  reference 
conditions  within  a regional  framework  (Davis  and  Simon  1995).  The  Montana  Water 
Quality  Division  recommends  the  collection  of  macroinvertebrate  samples  during  the 
summer  season  (June  21  to  September  21),  following  runoff.  This  is  because  the 
Montana  WQD  has  the  most  reference  data  from  this  season  and  because  stream  flows 
and  weather  conditions  are  most  suitable  for  field  work.  The  Montana  WQD  also 
recommends  that  for  monitoring  trends  at  a particular  site  over  time,  sampling  as  close 
as  possible  on  the  same  date  will  produce  the  best  results  (Bukantis  1995). 

Numerous  researchers  have  suggested  that  within-site,  within-stream. 


t>JCO{W 


uj«;.  ^ ' 't  ^^ib  * t*.-,.  €''*#  U S 'To  fV^w^04^ddl 

- ‘ J ‘ . <V  ■ ’ • 

-i»4Ci  •;,  is^w  .fi‘i.  1)  ^tu^ild  {go,  /A>to 


it:,'  .rjtf.qm  ’.’P  I-Air/''-  <•;  ■«»'■. Qqa©V  'i::'&?JiiSiff^Oi>'S^^  . jt 


'^J  ■- 

V-  2H,i  '•.  psL'  .;mmxs  '{Ua^^itr;  ■«  (4?-® 

‘ ,ar«i^3'3r:3  r^r.'nP  ^ 

' '4  • : ' f- 

« “S'' 

4>j.f:ti.v  ;v*  0*.  5‘^-viin^'TT  k-^^> atoPOfI 

» ! ' 

k’’  “ 

. n i/.-j  ij.Mss£i  -.'•*^‘rt8 

rrrjLvHad’^  v4^  . 

■ \ :'■ , ^ \ 

»-li  ^iru^ymKk75-iPj^o,4^il^^  -iiai^VflOWsJn^  ’*^S£w!3 


'<1  t: 


iL 


r»MT  f r |^4wa?joT  .(li  np4i^.i:$Cii!i'(ia 

’ • •*  - ' ’ . ' ' ' ' ' 

*'*  f ' ' ' 

/svoQ  3ic4,3iiB.'i^«aa5d  5(^t«  *’J'^)W  aar'^iiol.^ 

M ■'«1T  .:4fHiV';l«lf^  7r?l  ^y/?5l^2  8f»Gi^iUpou  i;«A. 

'.ifeXil  . iw  ^ J ■{.'»  ' ’liq  B !l  /«ii' a^rf«*s03>^; 


V *'** ' 


■A 


: .'s!  . diiv  'r.i  srfj  a.j^ti>iq  lUw  <*« 

_EH 
O™  , ^ 

,‘1!-  fc  ts  <ji^j 


15 


seasonal,  and  yearly  variability  of  stream  macroinvertebrate  communities  may  affect 
rapid  assessment  metrics  enough  to  make  them  unreliable  for  assessing  water  quality 
impairment  (Resh  and  Jackson  1993;  Bode  and  Novak  1995;  Angradi  and  Vinson 
1995).  To  be  effective,  rapid  assessment  metrics  must  be  able  to  discriminate  between 
impairment  and  namral  variability  (Bode  and  Novak  1995;  Barbour  et  al.  1995).  Very 
little  analyzed  data  exists  concerning  metric  variability  using  rapid  assessment  in 
mountain  streams  of  Montana. 

Macroinvertebrate  Variability  in  Mountain  Streams 

There  has  been  much  documentation  of  within-stream  and  seasonal  variability  of 
macroinvertebrates  in  mountain  streams  and  these  types  of  variability  are  well 
understood  concepts  (Ward  and  Kondratieff  1992;  Resh  and  Rosenberg  1989;  Yoder 
and  Rankin  1995).  For  example,  longitudinal  macroinvertebrate  community  changes 
within  a stream  from  headwaters  downstream  have  fostered  many  classification  and 
zonation  schemes  including  the  River  Continuum  Concept  (Vanote  et  al.  1980)  and 
lilies'  and  Botosaneanu's  (1963)  zonation  system. 

Seasonal  and  yearly  variability  studies  are  numerous.  Richards  et  al.  (1995) 
reported  a 3 to  4-  fold  increase  in  abundance  of  adult  Pteronarcys  californica 
(Plecoptera)  in  the  Madison  River  in  SW  Montana  between  1994  and  1995.  Gustafson 
(1990)  conducted  an  intensive  study  of  the  ecology  of  aquatic  insects  in  the  Gallatin 
River  drainage  in  SW  Montana.  From  his  work  on  the  mainstem  of  the  Gallatin  River, 
Gustafson  (1990)  reported  that  total  macroinvertebrate  density  was  very  low  the  first 
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summer  but  was  at  an  all  time  high  during  the  following  summer,  which  he  attributed 
to  an  early  spring  run-off  before  the  second  summer.  Gustafson  (1990)  also  reported 
that  in  autumn,  the  species  numbers  peak  due  to  the  appearance  of  species  with  winter 
life  cycles.  A winter  depression  of  species  numbers  occurred  due  to  winter  mortality 
or  downward  migration  from  the  top  layer  of  the  benthos.  The  peak  in  species  number 
during  the  spring  was  associated  with  the  appearance  of  species  with  summer  life 
cycles,  which  occurs  before  all  species  with  winter  cycles  have  emerged.  Gustafson 
(1990)  stated  that  spring  peaks  in  species  numbers  occurred  during  a depression  in  total 
invertebrate  density.  A decline  in  the  average  number  of  species  after  the  spring  peak 
was  associated  with  an  increased  total  density  because  many  of  the  species  with  winter 
life  cycles  were  absent  from  the  benthos  and  the  species  with  summer  life  cycles 
became  very  abundant  (Gustafson  1990).  Gustafson's  (1990)  conclusion  was  that  there 
were  low  macroinvertebrate  densities  in  summer,  then  an  increase  during  autumn,  a 
decline  in  winter,  a brief  increase  in  the  spring,  followed  by  a decline  to  the  summer 
low. 

The  State  of  Ohio  suggested  that  fundamental  decisions  need  to  be  made  early  in 
the  development  of  biocriteria,  including  accounting  for  seasonal  variability  (Yoder  and 
Rankin  1995).  The  State  of  Ohio  confines  their  collection  times  to  June  through 
September  and  is  presently  analyzing  temporal  variability  from  a special  Ohio  EPA 
methods  study  conducted  in  1981  (DeShon  1995).  Bode  and  Novak  (1995)  tested 
seasonal  variability  by  sampling  two  streams  monthly  for  a year.  They  suggested  that 
between-month  comparisons  should  not  be  done.  Reash  (1995)  reported  that  seasonal 
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and  yearly  variability  must  be  accounted  for  when  deriving  biocriteria.  Reash  (1995) 
based  this  by  monitoring  fish  populations  in  the  Ohio  River,  where  50%  of  the 
electrofishing  and  gill  net  samples  showed  significant  (5  % level)  temporal  variability  in 
catch  per  unit  effort  or  in  total  biomass.  Resh  and  Jackson  (1995)  reported  that  in  their 
study  of  the  accuracy  of  macroinvertebrate  rapid  assessment  measures,  16  of  the  22 
measures  tested  found  significant  (5%  level)  seasonal  differences.  A detailed 
knowledge  of  life  histories,  especially  voltinism  and  phenology,  of  macroinvertebrates 
encountered  in  biomonitoring  samples  should  help  alleviate  many  problems  associated 
with  seasonal  variability  (Rosenberg  and  Resh  1995;  Johnson  et  al.  1995). 

Functional  feeding  groups  are  also  known  to  follow  seasonal  patterns  of 
abundance  (Cummins  1974).  The  functional  feeding  group  classification  system 
"distinguishes  macroinvertebrate  taxa  that  perform  different  functions  within  aquatic 
ecosystems  with  respect  to  processing  of  nutritional  resource  categories"  (Merrit  and 
Cummins  1984).  Functional  feeding  groups  are  analogous  to  "guilds".  Many 
stoneflies  that  are  in  the  functional  feeding  group  'shredders'  have  short  life  cycles 
and  may  be  present  or  absent  from  samples  collected  at  certain  times  of  the  year 
(Gustafson  1990).  Gustafson  (1990)  reported  that  predators,  shredders,  and  collectors 
generally  increased  in  density  as  stream  order  increased,  while  the  density  of  scrapers 
peaked  in  the  mid-order  streams  in  the  Gallatin  River  drainage.  Scrapers  had  a high 
variance  at  all  stream  orders  and  headwater  streams  had  more  species  of  shredders  than 
streams  of  intermediate  size  (Gustafson  1990). 

At  present,  the  Montana  WQD  RBP's  for  the  mountain  stream  ecoregion  do 
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not  fully  address  within-site,  within-stream,  seasonal,  and  yearly  variability  of 
macroinvertebrate  communities,  although  variability  is  taken  into  account  by  setting 
wide  confidence  limits  on  impairment  (Bukantis  1996).  The  Montana  WQD  RBP's 
also  do  not  incorporate  the  potential  beneficial  uses  of  seasonal  variability  into  their 
assessments.  For  example,  the  higher  abundance  of  certain  taxa  at  different  seasons, 
other  than  summer,  that  may  be  more  sensitive  to  various  sources  of  impairment. 
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MATERIALS  AND  METHODS 
Within-site  Variability  Collection  Sites 

Twenty  six  sites  from  ten  mountain  streams  in  western  Montana  and  NW 
Wyoming  were  used  to  test  within-site  variability.  The  ten  mountain  streams  were 
Bozeman  Creek,  Hyalite  Creek,  Cottonwood  Creek,  Squaw  Creek,  Pebble  Creek,  Soda 
Butte  Creek,  North  Fork  Bear  Creek,  Middle  Fork  Bear  Creek,  Bear  Creek,  and 
Divide  Creek  (Figure  1).  A riffle  sequence  was  defined  as  a series  of  consecutive 
riffles  and  pools.  Three  macroinvertebrate  samples  were  collected  from  each  riffle 
sequence  (site),  one  sample  per  riffle,  in  1995.  Hyalite  Creek,  1994  and  Bozeman 
Creek,  August,  9,  1994  had  four  samples  collected  firom  a riffle  sequence  (1  sample 
per  riffle).  Collection  dates,  physical  and  chemical  data,  and  descriptions  of  the  streams 
are  presented  in  Table  2. 

Of  these  26  sites  ( riffle  sequence),  14  were  from  7 'least  impaired'  reference 
mountain  streams.  The  seven  reference  streams  included  5 sites  (riffle  sequences)  from 
Pebble  Creek,  (1995)  and  1 site  from  each  of  the  following  streams.  Cottonwood  Creek 
(1995),  N.  F.  Bear  Creek  (1995),  M.  F.  Bear  Creek  (1995),  Bear  Creek  (1995),  Soda 
Butte  Creek  (km  1,  1995),  Bozeman  Creek  (summer  1994)  and  Bozeman  Creek 
(summer  1995)(Table  2).  Mean,  variance,  standard  deviation  (s.d.),  coeffiecient  of 
variation  (CV)  and  95%  confidence  intervals  (Cl)  were  derived  for  each  metric  used. 
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Figure  1.  Map  of  study  sites 
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Table  2.  Descriptions  of  the  10  mountain  streams  sampled  in  this  study  including,  mountain  range,  drainage,  longitude, 

latitiude,  township,  range  1/4  section,  stream  order,  link  magnitude,  elevation  direction  of  stream  flow, 

pH,  conductivity,  Rosgen  stream  type,  estimated  discharge,  and  date  of  collection. 
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Table  2.  Continued. 
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Hyalite 

Creek 

B2c 

na 

June  27, 
1994 

Squaw 

Creek 

C4 

na 

July  12, 
1995 

Cotton- 

wood 

Creek 

C4 

na 

July  15, 
1995 

Soda 

Butte 

Creek 

B3,  B3 

na 

August 
16,  1995 

Divide 

Creek 

B3/B4 

na 

August 
30,  1995 

Bear 

Creek 

B3/B4 

0.35 

July  5, 
1995 

M.  F. 
Bear 
Creek 

B3c 

00 

July  5, 
1995 

N.  F. 
Bear 
Creek 

A2/A3b 

2.84 

July  5, 
1995 

Bozeman 

Creek 

C4b 

na 

Sept  20, 
1995 

Rosgen 

stream 



estimated 

discharge 

(m^'s) 

date 

sampled 

23 


Table  2.  continued.  Descriptions  of  Pebble  Creek  sites  sampled  in  this  study 
including,  mountain  range,  drainage,  longitude,  latitude,  township,  range  1/4 
section,  stream  order,  link  magnitude,  elevation  direction  of  stream  flow,  pH, 
conductivity,  Rosgen  stream  type,  estimated  discharge,  and  date  of  collection. 


measure 

Km  2 

Km  4 

Km  6.5 

Km  9 

Km  13 

Km  16 

Mountain 

range 

Absaroka 

range 

Absaroka 

range 

Absaroka 

range 

Absaroka 

range 

Absaroka 

range 

Absaroka 

range 

drainage 

Lamar 

River 

(Yellow- 

stone 

River) 

Lamar 

River 

Lamar 

River 

Lamar 

River 

Lamar 

River 

Lamar 

River 

longitude 

110"  02’ 

na 

na 

na 

na 

110"  10' 

latitude 

45"  or 

na 

na' 

na 

na 

44"  13’ 

township 
range 
1/4  section 

T9S/R14E 
sec  19 

na 

na 

na 

na 

T57N/R1 1 
OW 
sec  13 

stream 

order 

2 

2 

3 

3 

3 

3 

link 

magnitude 

5 

6 

8 

10 

17 

23 

elevation 

(meters) 

2733 

2667 

2600 

2533 

2483 

2300 

direction 
of  stream 
flow 

SW 

SW 

S 

S 

S 

S 

pH 

8.8 

8.6 

8.6 

8.0 

8.2 

8.5 

conduct- 

ivity 

(micro 

siemens) 

170 

190 

200 

170 

170 

160 

Rosgen 

stream 

type 

B3c 

B3c 

C3/C4 

C3/C4 

B2/B3 

B2/B3 

estimated 

discharge 

(mVs) 

1.08 

na 

2.64 

2.73 

na 

na 

date 

sampled 

Augtust  16, 
1995 

August  16, 
1996 

August  15, 
1995 

August  15, 
1995 

August  15, 
1995 

August  16, 
1995 
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Three  samples  collected  from  N.  F.  Bear  Creek,  M.  F.  Bear  Creek,  and  Bear 
Creek  (collectively  referred  to  as  the  Three  Bear  Creeks  in  this  study)  were  collected 
on  July  5,  1995,  while  one  sample  was  collected  from  the  same  riffle  on  June  24,  1994 
from  each  stream.  The  three  Bear  Creeks  originate  in  the  Madison  range  of  SW 
Montana  within  the  Lee  Metcalf  Wilderness  (Table  2 and  Figure  1).  Sample  sites  were 
within  the  wilderness  area  and  were  assumed  to  be  'least  impaired'  sites. 

Divide  Creek  is  a third  order  stream  near  St.  Mary,  MT  and  forms  part  of  the 
eastern  border  between  Glacier  National  Park  and  the  Blackfeet  Indian  Reservation  in 
NW  Montana  (Figure  1.  and  Table  2 ).  Divide  Creek  has  been  drastically  channelized 
in  an  attempt  to  prevent  flooding  of  private  businesses  in  St.  Mary  and  to  prevent  a 
historic  bridge  (operated  by  the  National  Park  Service)  from  being  washed  out  (Rosgen 
1994b).  Divide  Creek  is  usually  channelized  by  bulldozers  yearly  (Rosgen  1994b). 
Therefore,  Divide  Creek  was  used  as  an  'impaired'  site. 

One  sample  was  also  collected  on  October  20,  1994  from  Squaw  Creek.  Squaw 
Creek  has  had  extensive  logging  associated  activities  in  its  drainage  including  the 
removal  of  most  of  the  large  woody  debris  in  the  lower  mainstem  where  the  samples 
were  collected  (Table  3)  and  was  considered  'impaired'. 

One  sample  from  two  sites  (km  1 and  km  2)  was  also  collected  on  August  13, 
1994  from  Soda  Butte  Creek.  Soda  Butte  Creek  is  located  near  Cooke  City,  SW 
Montana  (Table  2 and  Figure  1).  Soda  Butte  Creek  receives  mine  leachate  from  the 
McLaren  tailings  in  Cooke  City  before  it  flows  into  Yellowstone  National  Park.  One 
collection  site  was  located  above  the  McLaren  tailings  (km  1)  and  considered  'least 
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impaired'  and  one  site  was  located  directly  below  the  McLaren  tailings  (km  2)  and  was 
obviously  'impaired'.  The  amount  of  logging  associated  activity  in  Bozeman  Creek, 
Hyalite  Creek,  Squaw  Creek,  and  Cottonwood  Creek  are  in  Table  3. 


Table  3.  Number  of  road  kilometers  per  square  kilometer  (km/km^),  number  of  clear- 
cut  hectares  per  square  kilometer  (hect/km^),  and  the  number  of  cubic  meters  of  woody 
debris  per  stream  kilometer  (m^/km)  in  the  Bozeman,  Hyalite,  Cottonwood,  and  Squaw 
Creek  drainages^ 


Bozeman 

Creek 

Hyalite  Creek 

Cottonwood 

Creek 

Squaw  Creek 

roads 

(km/km^) 

0.86 

1.15 

0.24 

1.26 

clear-cuts® 
(hect./  km^) 

7.71 

12.52 

4.96 

23.28 

woody  debris‘s 
(mVkm) 

152.01 

66.52 

191.03 

46.18 

® data  was  not  available  on  the  age  of  the  clear-cuts  or  how  close  to  the  stream  the 
dear-cuts  occured. 

**  woody  debris  was  measured  along  a 3.2  km  (2  mile)  similar  section  of  each  stream. 


Within-stream  Variability  Collection  Sites 


Six  collection  sites  (3  samples  per  site)  in  Pebble  Creek  (Y ellowstone  National 
Park,  Wyoming  and  Montana)  were  used  to  test  within-stream  variability  (Figure  1 and 
Table  2).  These  sites  were  also  used  to  test  within-site  variability.  Pebble  Creek  was 
considered  to  be  a 'least  impaired'  mountain  stream.  Any  changes  in  the  metrics 
calculated  for  Pebble  Creek  were  assumed  due  to  natural  causes.  Collection  sites  in 


Pebble  Creek  were  at  approximately  2 km,  4 km,  6.5  km,  9 km,  13  km,  and  16  km 
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from  its  headwaters.  Pebble  Creek  sites  were  labelled  site  2 km,  site  4 km,  site  6.5  km 
and  so  forth  for  future  reference.  Principal  Components  Analysis  using  raw  species 
abundances  from  each  sample  (3  per  site)  was  used  to  explore  longitudinal  relationships 
between  each  site.  Mean,  variance,  standard  deviation  (s.d.),  coefficient  of  variation 
(CV),  and  95%  confidence  intervals  (Cl)  for  twenty  metrics  were  calculated  for  each 
site. 

Six  collections  sites  in  Soda  Butte  Creek  (Montana  and  Yellowstone  National 
Park,  Wyoming)  were  also  used  to  test  within-stream  variabililty  (Figure  1)  using  PC  A 
of  raw  species  abundance  data  for  each  sample  as  a sampling  unit.  Soda  Butte  Creek  is 
an  'impaired'  stream,  therefore,  within-stream  variability  analysis  was  complicated  by 
its  impairment.  Sites  were  labelled  for  their  distance  from  Soda  Butte  Creek's  origin  at 
1km,  2 km,  4 km,  8 km,  11  km  and  16  km. 
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Seasonal  Variability  Collection  Sites 

Samples  used  to  test  for  seasonal  variability  were  collected  from  the  same  riffle 
sequence  (site)  from  Bozeman  Creek  (Gallatin  River  drainage,  SW  Montana)  at  five 
different  times;  four  samples  on  August  9,  1994,  three  samples  on  October  21,  1994, 
three  samples  from  February  5,  1995,  three  samples  on  May  5,  1995,  and  three 
samples  on  September  20,  1995. 

Bozeman  Creek,  also  known  as  Sourdough  Creek,  near  Bozeman,  MT  (Figure  1 
and  Table  2)  is  part  of  the  city  of  Bozeman's  water  supply.  There  is  an  unimproved 
road  along  Bozeman  Creek  for  much  of  its  length  which  is  used  by  recreationalists. 
Bozeman  Creek  is  closed  to  unauthorized  vehicles.  Only  limited  logging  activity  has 
occured  in  the  Bozeman  Creek  drainage  (Table  3).  The  stream  corridor  recieves  limited 
use  and  is  protected  from  logging  associated  activities.  Because  it  is  a city  water 
supply,  Bozeman  Creek  was  used  as  a reference  or  'least  impaired'  stream. 

Macroinvertebrate  Collection 

Macroinvertebrate  samples  were  collected,  sorted,  and  taxonomically  identified 
using  the  Water  Quality  Division  (WQD),  State  of  Montana,  Standard  Operating 
Procedures,  except  that  no  sub-sampling  was  done  (State  of  Montana  recommends  a 
300  organism  sub-sample).  A one-minute  traveling  kick  using  a standard  D-frame  net 
with  a 1 mm  mesh  size  was  used  in  riffle  habitats  for  all  samples  collected  (Bukantis 
1995).  Riffle  samples  were  collected  by  vigorously  kicking  the  substrate  and  moving 
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diagonally  downstream  and  across  the  riffle  in  an  effort  to  sample  all  of  the  riffle  sub- 
habitats at  an  approximate  rate  of  0.3  meters  (1  foot)  per  second.  Timing  of  the  sample 
collection  for  one  minute  provided  for  a semi-quantitative  analysis  used  in  the 
computation  of  various  metrics. 

Samples  were  then  either  placed  in  a shallow  pan  and  all  macroinvertebrates 

hand-picked  and  placed  in  95  % ethanol  at  the  collection  site  or  all  of  the  contents 

collected  in  the  net,  including  debris,  were  directly  placed  in  95%  ethanol  and  sorted  at 

a later  time.  For  the  purpose  of  this  study  all  of  the  organisms  collected  in  the  samples 

were  sorted  and  identified  regardless  of  their  numbers.  Macroinvertebrates  were 

identified  to  the  taxonomic  level  suggested  by  the  State  of  Montana,  usually  genus  or 

species,  by  using  taxonomic  keys  from  Merrit  and  Cummins  1983,  Stewart  and  Stark 

1993,  Pennak  1989,  Wiggins  1977,  and  Gustafson  1994.  More  difficult  identifications 

* 

were  verified  by  Dr.  Dan  Gustafson  at  Montana  State  University.  Voucher  specimens 
were  placed  in  the  Montana  Entomology  Collection  at  Montana  State  University. 

Physical  and  Chemical  Data 

Physical  and  chemical  stream  data  were  collected  at  the  time  of  the 
macroinvertebrate  collections  at  each  site  and  included:  1)  temperature,  2)  pH,  3) 
conductivity,  4)  gradient,  5)  velocity,  6)  width,  7)  depth,  and  8)  substrate  composition. 
Water  temperature  was  measured  with  a calibrated  mercury  bulb  thermometer.  pH  was 
measured  with  a Cole-Parmer  pH  Tester  2™.  Conductivity  was  measured  with  a Cole- 
Parmer  TDS  Tester™.  Stream  gradient  at  the  riffle  site  was  estimated  using  a 60  meter 
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line  and  line  level.  Velocity  was  estimated  by  floating  a stick  for  a distance  of  60 
meters  centered  on  the  riffle  section  and  timed  on  a watch  to  the  nearest  second.  Five 
passes  were  made  and  then  averaged.  A correction  factor  of  0.85  was  used  to  adjust 
for  average  velocity  within  the  water  column  (Harrelson  et  al.  1994).  Stream  width  was 
measured  at  estimated  bank  full  stage  to  the  nearest  tenth  of  a meter  and  stream  depth 
was  measured  at  estimated  bank  full  stage  to  the  nearest  centimeter  at  ten  intervals  and 
then  averaged.  Velocity,  width,  and  average  depth  were  used  to  roughly  estimate 
discharge  at  sites  where  no  gaging  stations  are  established.  Substrate  composition  was 
measured  using  a Wolman  100  pebble  count  (Wolman  1954;  Kondolf  and  Li  1992). 
Measurements  of  gradient  and  substrate  size  composition  along  with  an  estimate  of 
sinuosity  and  entrenchment  were  used  to  classify  each  study  site  using  the  Rosgen 
stream  classification  system  (Rosgen  1994a). 

Stream  order,  link  magnitude,  elevation,  longitude,  latitude,  township,  range, 
1/4  section,  mountain  range,  major  river  drainage,  and  flow  direction  were  made  from 
U.S.G.S.  15-  minute  topographic  maps  and  U.S.D.A.  Forest  Service  maps.  Stream 
order  was  classified  using  the  following  method;  a first  order  stream  is  the  headwater 
stream,  a second  order  stream  is  formed  where  two  first  order  streams  join,  a third 
order  stream  is  formed  when  two  second  order  streams  join  and  a fourth  order  stream 
is  formed  when  two  third  order  streams  join  (Strahler  1957).  Link  magnitude  is  the 
total  number  of  first  order  streams  entering  within  the  drainage  upstream  of  the  study 
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Metrics  and  Statistical  Analysis 

Metrics  calculated  and  used  in  this  study  are  listed  in  Table  1 . Mean,  95  % 
confidence  interval  (Cl),  variance,  standard  deviation  (s.d.)  and  coefficient  of  variation 
(CV)  were  calculated  for  sites  with  more  than  one  sample  per  site  using  QuatroPro 
(Borland  1993).  All  (+. ) intervals  reported  in  this  study  are  95%  confidence  intervals 
(Cl).  Richness,  evenness,  and  diversity  indices  were  computed  for  most  samples  using 
Statistical  Ecology  (Ludwig  and  Reynolds  1988).  Principal  Components  Analysis 
(PCA)  was  used  to  explore  variability  in  the  macroinvertebrate  community  in  Bozeman 
Creek,  Pebble  Creek,  Soda  Butte  Creek,  and  the  three  Bear  Creeks  by  using  the  raw 
species  abundance  data  for  each  sample  as  a sampling  unit  using  Statistical  Ecology 
(Ludwig  and  Reynolds  1988)  and  QuatroPro  (Borland  1993),  Typically  the  first  three 
eigenvectors  describe  most  of  the  variability  (Ludwig  and  Reynolds  1988),  therefore, 
only  the  first  three  eigenvector  values  were  computed. 
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RESULTS 

Within-site  Variability 

Mean  within-site  CV's  ranged  from  a low  of  0.05  for  % EPT  taxa  to  a high  of 
0.43  for  % Chironomidae  (Table  4).  The  second  highest  mean  CV  value  was  the 
number  of  organisms  metric  at  0.23.  The  other  19  metrics  CV  values  were  less  than 
0.23. 

Longitudinal  comparisons  of  mean  and  95%  Cl's  for  15  metrics  from  Pebble 
Creek  also  showing  within-site  variability  are  in  the  within-stream  portion  of  the 
Results  section.  Seasonal  mean  and  95%  Cl's  showing  within-site  variability  for  17 
metrics  from  Bozeman  Creek  are  in  the  seasonal  portion  of  the  Results  section. 

Voltinism  metrics  were  not  used  for  further  analysis  because  the  range  in 
variability  calculated  for  within-sites  in  Bozeman  Creek  was  considered  too  large 
(range  of  0.00  CV  to  1.73  CV)  and  also  because  present  knowledge  of  voltinism  in 


mountain  stream  macroinvertebrates  in  Montana  is  limited. 
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Table  4.  Mean  of  CV's  for  reference  stream  sites  (n  = 34  samples,  except  Metals 
Tolerance  Index:  n = 17  samples). 


metric 

mean 

CV 

it  of  organisms 

0.23 

taxa  richness’ 

0.09 

EPT  richness^ 

0.08 

it  of  Plecoptera  taxa 

0.18 

it  of  Ephemeroptera  taxa 

0.07 

it  of  Trichoptera  taxa 

0.12 

a of  shredder  taxa 

0.12 

it  of  scraper  taxa 

0.09 

a of  eg  and  fg  taxa^ 

0.09 

it  of  predator  taxa 

0.22 

% shredder  taxa 

0.17 

% scraper  taxa 

0.11 

%cg  and  fg  taxa 

,0.11 

% predator  taxa 

0.19 

% dominant 

0.18 

HBI' 

0.15 

Shannon  Diversity  Index 

0.07 

%EPT  richness 

0.05 

%EPT  organisms 

0.06 

% Chironomidae 

0.43 

Metals  Tolerance  Index 

0.15 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 


Within-stream  Variability 

Pebble  Creek 

Principal  Components  Analysis  of  the  6 sites  from  Pebble  Creek  using  raw 
species  abundance  data  showed  that  macroinvertebrate  samples  from  km  2 (N  = 3),  the 
most  upstream  headwater  site  (2nd  order  stream  segment),  clustered  together  and  were 
distinct  from  the  other  sites.  Site  km  4 (N  = 1)  also  a 2nd  order  stream  segment, 
clustered  separately  from  the  rest  of  the  sites.  Macroinvertebrate  samples  from  sites 
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km  6.5,  km  9,  km  13,  and  km  16  (N  = 3)  formed  a loose  cluster  with  sites  km  13  and 
km  16  clustering  more  tightly  than  the  other  sites  (Table  5 and  Figure  2). 


Table  5.  Eigenvectors  defining  the  first  three  principal  components  from  16  samples 
from  Pebble  Creek,  1995. 


Sample  Unit I 


Km  2 

1 

-0.16 

2 

0.07 

3 

-0.07 

Km  4 

4 

-0.7 

Km  6.5 

5 

-0.11 

6 

0.57 

7 

-0.66 

Km9 

8 

-0.82 

9 

-0.52 

10 

2.69 

km  13 

11 

0 

12 

-0.23 

13 

-0.11 

km  i6 

14 

0.52 

15 

0.04 

16 

0.11 

II 

-1.22 

-0.95 

-1.14 

■-0"67 

0"54 

1.29 

1.29 

oil' 

0.5 

0.02 

••••••••••••I 

0.14 

0.05 

0.05 

■0.'06 

-0.21 

-0.03 


III 

0.37 

-0.3 

0.05 

■‘i'.'68 

■o."i'6 

0.34 

0.57 

'0.'08 

0.05 


-0.63 

-0.48 

-0.'56 

-0.85 

-0.75 


The  mean,  variance,  standard  deviation.  Cl,  and  CV  values  for  the  six  Pebble 


Creek  sites  are  in  Tables  6 through  10. 
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Figure  2.  Principal  Components  Analysis  for  Pebble  Creek. 
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Table  6.  Mean  (+.  95%  Cl),  variance,  standard  deviation  (s.  d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Pebble  Creek  1995  site  km  2 (N=3). 


metric 

mean 

(±95% 

C.I.) 

variance 

s.d. 

CV 

# of  organisms 

588.33 
(±  243.73) 

46386.00 

215.37 

0.38 

taxa  richness' 

25.33 
(±  3.46) 

9.33 

3.06 

0.12 

EPT  richness^ 

17.67 

(±2.85) 

6.33 

2.52 

0.14 

# of  Plecoptera  taxa 

3.33 

(±1.30) 

1.33 

1.15 

0.35 

# of  Ephemeroptera 
taxa 

7.67 

(±1.30) 

1.33 

1.15 

0.15 

M of  Trichoptera  taxa 

6.67 

(±0.66) 

0.33 

0.58 

0.09 

# of  shredder  taxa 

1.00 

(±0.00) 

0.00 

0.00 

0.00 

if  of  scraper  taxa 

3.33 

(±0.65) 

0.33 

0.58 

0.17 

it  of  eg  and  fg  taxa-' 

5.67 

(+1.31) 

1.33 

1.15 

0.20 

a of  predator  taxa 

7.67 

(±  1.73) 

2.33 

1.53 

0.20 

% shredder  taxa 

4.00 

(±1.00) 

0.00 

0.00 

0.12 

% scraper  taxa 

13.00 

(±2.00) 

0.00 

1.00 

0.11 

% eg  + fg  taxa 

21.00 

(±3.00) 

0.00 

3.00 

0.03 

% predator  taxa 

29.00 

(±4.00) 

0.00 

4.00 

0.13 

% dominant 

26.38 

(±11.08) 

95.95 

9.80 

0.37 

HBP 

4.89 

(±0.89) 

0.62 

0.79 

0.16 

Shannon  Diversity 
Index 

2.48 

(±0.19) 

0.03 

0.17 

0.07 

%EPT  richness 

67.11 

na 

na 

na 

MTI 

2.20 

(±  0.25) 

0.05 

0.22 

0.10 

EPT  taxa 

17.67 

(±2.85) 

6.33 

2.52 

0.14 

%EPT  organisms 

39.00 
(±  13.00) 

1.00 

11.00 

0.28 

% Chironomidae 

54.00 
(±  14.00) 

2.00 

12.00 

0.23 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeropiera,  Plecoptera,  and  Trichoptera  taxa  in  a sample. 

^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa;  * HBI  is  Hilsenhoff  Biotic  Index 
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Table  7.  Mean  (+.95%  C.I.)^  variance,  standard  deviation  (s.  d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Pebble  Creek,  1995  site  km  6.5  (N  = 3). 


metric 

mean 

(±95% 

C.I.) 

variance 

s.  d. 

CV 

# of  organisms 

429.00 
C±  137.49) 

14763.00 

121.50 

0.28 

taxa  richness’ 

28.33 
(±  2.85) 

6.33 

2.52 

0.09 

EPT  richness^ 

21.00 

(+1.13) 

1.00 

1.00 

0.05 

# of  Plecoptera 
taxa 

4.67 

(+0.66) 

0.33 

0.58 

0.12 

#of 

Ephemeroptera 

taxa 

10.33 
(±  0.66) 

0.33 

0.58 

0.06 

# of  Trichoptera 
taxa 

6.00 

(+0.00) 

0.00 

0.00 

0.00 

# of  shredder 
taxa 

1.00 

(+0.00) 

0.00 

0.00 

0.00 

# of  scraper  taxa 

5.67 

(+0.65) 

0.33 

0.58 

0.10 

# of  eg  and  fg 
taxa^ 

7.67 

(±1.73) 

2.33 

1.53 

0.20 

# of  predator 
taxa 

7.67 

(+1.31) 

1.33 

1.15 

0.15 

% shredder  taxa 

4.00 

(±0.00) 

0.00 

0.00 

0.09 

% scraper  taxa 

20.00 

(±3.00) 

0.00 

3.00 

0.13 

%cg  and  fg  taxa 

28.00 

(±9.00) 

1.00 

8.00 

0.28 

% predator  taxa 

27.00 

(±5.00) 

0.00 

5.00 

0.18 

% dominant 
taxa 

13.50 
(±  1.47) 

1.69 

1.30 

0.10 

0.78 

(±0.03) 

0.00 

0.03 

0.04 

Shannon 
Diversity  Index 

2.74 

(±  0.02) 

0.00 

0.02 

0.01 

%EPT  richness 

74.30 
(±  3.43) 

0.09 

3.03 

0.04 

% EPT 
organisms 

93.00 
(±  3.00) 

0.00 

3.00 

0.03 

% 

Chironomidae 

4.00 

(±2.00) 

0.00 

2.00 

0.44 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  8.  Mean  (i  95%  Cl),  variance,  standard  deviation  (s.  d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Pebble  Creek,  site  km  9 (N  = 3). 


metric 

mean 

(±  95%  C.I.) 

variance 

s.d. 

CV 

a of  organisms 

483.67 
C±  63.26) 

3124.33 

55.90 

0.12 

taxa  richness' 

28.33 

C+3.64) 

10.33 

3.21 

0.11 

EPT  richness^ 

22.67 

(+1.31) 

1.33 

1.15 

0.05 

it  of  Plecoptera 
taxa 

5.00 

(+1.13) 

1.00 

1.00 

0.20 

#of 

Ephemeroptera 

taxa 

11.00 

(±1.13) 

1.00 

1.00 

0.09 

it  of  Trichoptera 
taxa 

6.67 

(+1.73) 

2.33 

1.53 

0.23 

a of  shredder  taxa 

1.33 

(+0.65) 

0.33 

0.58 

0.43 

a of  scraper  taxa 

5.00 

(+0.00) 

0.00 

0.00 

0.00 

a of  eg  and  fg 
taxa^ 

8.33 

(+0.65) 

0.33 

0.58 

0.07 

it  of  predator  taxa 

7.33 

(+1.31) 

1.33 

1.15 

0.16 

% shredder  taxa 

5.00 

(+2.00) 

0.00 

0.01 

0.31 

% scraper  taxa 

18.00 

(+2.00) 

0.00 

2.00 

0.11 

% eg  and  fg  taxa 

30.00 

(+5.00) 

0.00 

4.00 

0.14 

% predator  taxa 

26.00 

(+4.00) 

0.00 

3.00 

0.13 

% dominant  taxa 

27.40 

(+6.89) 

37.03 

6.08 

0.22 

iffip 

1.65 

(±0.06) 

0.00 

0.05 

0.15 

Shannon  Diversity 
Index 

2.56 
(±  0.09) 

0.01 

0.08 

0.03 

%EPT  richness 

80.37 
(±  5.55) 

0.24 

4.90 

0.06 

%EPT  organisms 

89.00 
(±  2.00) 

0.00 

2.00 

0.02 

% Chironimidae 

9.00 

(±  2.00) 

0.00 

2.00 

0.22 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  9.  Mean  (±95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Pebble  Creek,  site  km  13  (N=3). 


metric 

mean 

(±  95%  C.I.) 

variance 

s.d. 

CV 

a of  organisms 

241.67 

(±39.14) 

1196.33 

34.59 

0.14 

taxa  richness' 

23.67 

(±2.61) 

5.33 

2.31 

0.10 

EPT  richness^ 

18.33 

(±1.31) 

1.33 

1.15 

0.06 

# of  Plecoptera 
taxa 

2.33 

(±  0.66) 

0.33 

0.58 

0.25 

#of 

Ephemeroptera 

taxa 

9.67 

(±1.30) 

1.33 

1.15 

0.12 

if  of  Trichoptera 
taxa 

6.33 

(±0.66) 

0.33 

0.58 

0.09 

if  of  shredder  taxa 

1.00 

(±0.00) 

0.00 

0.00 

0.00 

a of  scraper  taxa 

5.33 

(±0.65) 

0.33 

0.58 

0.11 

it  of  eg  and  fg 
taxa^ 

6.67 

(±0.65) 

0.33 

0.58 

0.09 

it  of  predator  taxa 

5.00 

(+1.13) 

1.00 

1.00 

0.20 

% shredder  taxa 

4.00 

(±0.00 

0.00 

0.00 

0.00 

% scraper  taxa 

23.00 

(±3.00) 

0.00 

2.00 

0.10 

% eg  + fg  taxa 

28.00 

(±0.00) 

0.00 

0.00 

0.01 

% predator  taxa 

21.00 

(±5.00) 

0.00 

5.00 

0.21 

% dominant  taxa 

18.73 

(±2.36) 

4.36 

2.09 

0.11 

HBP 

1.04 

(±0.28) 

0.06 

0.25 

0.24 

Shannon  Diversity 
Index 

2.61 

(+0.10) 

0.01 

0.09 

0.03 

% EPT  taxa 

77.44 

(+6.10) 

0.25 

4.81 

0.08 

%EPT  organisms 

89.00 

(+3.00) 

0.00 

3.00 

0.03 

% Chironimidae 

5.00 
(+  3.00) 

0.00 

3.00 

0.54 

’Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  10.  Mean  (+.95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Pebble  Creek,  site  km  16  (N=3). 


metric 

mean 

(±95%  C.I.) 

variance 

s.d. 

CV 

it  of  organisms 

156.33 

(±51.32) 

2056.33 

45.35 

0.29 

taxa  richness' 

24.67 

(±1.73) 

2.33 

1.53 

0.06 

EPT  richness' 

21.67 

(+1.73) 

2.33 

1.53 

0.07 

a of  Plecoptera 
taxa 

4.67 

(±0.66) 

0.33 

0.58 

0.12 

#of 

Ephemeroptera 

taxa 

11.00 

(±1.13) 

1.00 

1.00 

0.09 

it  of  Trichoptera 
taxa 

6.00 

(±1.13) 

1.00 

1.00 

0.17 

if  of  shredder  taxa 

1.00 

(±0.00) 

0.00 

0.00 

0.00 

it  of  scraper  taxa 

3.67 

(±  0.65) 

0.33 

0.58 

0.16 

it  of  eg  and  fg 
taxa^ 

9.67 

(±  0.65) 

0.33 

0.58 

0.06 

it  of  predator  taxa 

7.00 

(±0.00) 

0.00 

0.00 

0.00 

% shredder  taxa 

4.00 

(±0.00) 

0.00 

0.00 

0.06 

% scraper  taxa 

15.00 

(±2.00) 

0.00 

2.00 

0.11 

%cg  and  fg  taxa 

39.00 

(±4.00) 

0.00 

4.00 

0.10 

% predator  taxa 

28.00 

(±2.00) 

0.00 

2.00 

0.06 

% dominant  taxa 

11.70 

(±1.46) 

1.66 

1.29 

0.11 

ifflP 

1.11 

(±0.42) 

14.00 

37.00 

0.33 

Shannon  Diversity 
Index 

2.89 

(±  0.08) 

0.00 

0.07 

0.02 

% EPT  richness 

87.84 
(±  3.67) 

1.39 

2.43 

0.08 

%EPT  organisms 

82.00 

(±3.00) 

0.00 

2.00 

0.03 

% Chironirnidae 

11.00 
(+  2.00) 

0.00 

0.02 

0.15 

‘Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Using  the  mean  and  95  % Cl  for  within-stream  metric  comparison,  the  number 
of  organisms  per  sample  in  Pebble  Creek  decreased  from  headwaters  downstream,  with 
km  2 being  most  variable  (Figure  3).  Taxa  richness  was  highest  in  the  mid-reach  C- 
type  channels  (Rosgen  method),  km  6.5,  and  km  9 (Figure  4).  Number  of  EPT  (Figure 
5)  taxa  increased  downstream.  Number  of  Ephemeroptera  taxa  increased  in  a 
downstream  direction  (Figure  6).  Number  of  Plecoptera  taxa  was  highest  in  the  mid- 
reach C-type  channels,  km  6.5  and  km  9 (Figure  7).  Number  of  Trichoptera  taxa  was 
fairly  constant  throughout  Pebble  Creek,  with  km  4 having  the  highest  value  of  9 taxa 
(Figure  8).  The  number  of  shredder  taxa  was  consistent  at  one  taxon  with  little 
variability  (Figure  9).  Number  of  scraper  taxa  was  highest  in  the  mid-reaches  of  Pebble 
Creek  (sites  km  6.5  and  km  9)(Figure  10).  Number  of  collector-gatherer  and  fiter- 
gatherer  taxa  increased  from  the  upstream  to  downstream  sites  (Figure  11).  Number  of 
predator  taxa  did  not  show  any  trend  for  the  sites  (Figure  12).  Percent  shredder  taxa 
was  3%  to  7%  (95%  Cl)  throughout  Pebble  Creek  (Figure  13). Percent  scraper  taxa 
increased  downstream  except  for  a decrease  at  km  16  (Figure  14).  Percent  collector- 
gatherer  and  filterer-gatherer  taxa  increased  downstream  (Figure  15).  Percent  predator 
taxa  was  fairly  constant  throughout  (Figure  16).  Percent  dominant  taxa  was  highly 
variable  within-sites  and  between  sites  (Figure  17).  Percent  Chironomidae  organisms 
was  highest  and  most  variable  at  km  2 and  decresed  markedly  downstream  (Figure  18). 
Percent  EPT  taxa  increased  sharply  from  km  2 downstream  (Figure  19).  Shannon 
Diversity  index  also  tended  to  increase  downstream  (Figure  20).  HBI  was  highest  at 
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km  2 and  decreased  rapidly  downstream  with  little  within-site  variability  (Figure  21). 

The  large  values  of  HBI  and  MTI  in  site  Km  2 and  Km  4 were  due  to  the  high 
relative  abundance  of  Chironomidae  larvae  that  were  identified  to  genus  and  sub-family 
levels.  HBI  values  for  Chironomidae  larvae  at  the  genus  and  sub-family  level  are 
based  on  data  from  pollution  tolerant  species  and  may  not  reflect  pollution  tolerances  of 
mountain  stream  species  of  Chironomidae  larvae  within  these  genera  and  sub-families. 
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Figure  3.  Mean  (±95%  Cl)  number  of  organisms  for  six  sites  from  Pebble  Creek. 
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Figure  4.  Mean  (±95%  Cl)  number  of  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  5.  Mean  (±95%  Cl)  number  of  EPT  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  6.  Mean  C±  95%  Cl)  number  of  Ephemeroptera  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  7.  Mean  (±95%  Cl)  number  of  Plecoptera  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  8.  Mean  C±.  95%  Cl)  number  of  Trichoptera  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  10.  Mean  (±95%  Cl)  number  of  scraper  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  1 1 Mean  C±.  95  % Cl)  number  of  collector-gatherer  and  filter-gatherer  taxa  for  six  sites  from  Pebble  Creek. 


50 


10  %96  PUB  UB9LU 


BXBl  pUB  60  jO  # 


km  2 km  4 km6.5  km9  km13  km16 


I 


1 


ti 


ii-l.'" 


Figure  12.  Mean  C±  95%  Cl)  number  of  predator  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  13.  Mean  (±95%  Cl)  % shredder  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  14.  Mean  (±.95%  Cl)  % scraper  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  15.  Mean  (±95%  Cl)  % collector-gatherer  and  filter-gatherer  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  16.  Mean  (±95%  Cl)  % predator  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  18.  Mean  C±  95%  Cl)  % Chironomidae  organisms  for  six  sites  from  Pebble  Creek. 
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Figure  19.  Mean  C±  95%  Cl)  % EPT  taxa  for  six  sites  from  Pebble  Creek. 
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Figure  20.  Mean  (±95%  Cl)  Shannon  Diversity  Index  values  for  six  sites  from  Pebble  Creek. 
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Figure  21.  Mean  (i  95%  Cl)  HBI  values  for  six  sites  from  Pebble  Creek. 
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Seasonal  Variability 

Bozeman  Creek 

Results  of  Principal  Components  Analysis  using  raw  species  abundance  data 
from  each  sample  as  a sampling  unit  from  each  of  four  seasons  for  Bozeman  Creek  are 
in  Table  11.  Graphical  analysis  using  the  first  two  principal  component  eigenvalues  as 
axis'  show  that  Winter  1995  and  Spring  1995  were  most  similar  (Figure  22). 


Table  11.  Eigenvectors  defining  the  first  three  principal  components  from  thirteen  riffle 
samples  from  Bozeman  Creek  for  four  seasons  using  raw  species  abundance  data  from 
each  sample  as  a sampling  unit. 


Sample  Unit 

I 

n 

m 

Summer 

1 

-1.451 

-1.05 

-0.558 

2 

-0.714 

-0.399 

-0.013 

3 

-1.368 

-0.814 

-0.85 

4 

-0.986 

-0.877 

-0.321 

Autumn 

5 

-0.7 

1.988 

0.522 

6 

-0.325 

0.924 

0.091 

7 

-0.155 

1.776 

-0.061 

Winter 

8 

0.855 

-0.148 

-0.047 

9 

0.76 

-0.04 

-1.7 

10 

0.523 

-0.08 

-0.725 

Spring 

11 

0.595 

-0.445 

-0.523 

12 

1.773 

-0.614 

1.501 

13 

1.193 

-0.222 

0.342 

The  mean  (+.95%  C.I.),  variance,  standard  deviation  (s.  d.)  and  coefficient  of 
variation  (CV)  for  21  metrics  calculated  for  the  different  season  samples  from  Bozeman 
Creek  are  presented  in  Tables  12  through  16. 
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Figure  22.  Principal  components  analysis  for  Bozeman  Creek 
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Table  12.  Mean  (+.  95%  Cl),  variance,  standard  deviation  (s.  d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Bozeman  Creek,  August  9,  1994  (N=4). 


metric 

mean 

(+  95%  C.I.) 

variance 

s.  d. 

CV 

if  of  organisms 

315.75 
(±  63.56) 

4206.90 

64.86 

0.21 

taxa  richness' 

27.75 

(+1.24) 

1.58 

1.26 

0.05 

EPT  richness^ 

19.75 

(+1.47) 

2.25 

1.50 

0.08 

it  of  Plecoptera 
taxa 

2.00 

(±0.00) 

0.00 

0.00 

0.00 

^^'of 

Ephemeroptera 

taxa 

9.25 

(±0.94) 

0.92 

0.96 

0.10 

it  of  Trichoptera 
taxa 

8.5 

(±0.57) 

0.33 

0.58 

0.07 

it  of  shredder  taxa 

2.25 

(±0.96) 

0.92 

0.98 

0.44 

it  of  scraper  taxa 

4.00 

(±0.80) 

0.67 

0.82 

0.21 

it  of  eg  and  fg 
taxa 

10.5 

(±0.57) 

0.33 

0.58 

0.06 

it  of  predator  taxa 

5.00 

(±1.38) 

2.00 

1.41 

0.28 

% shredder  taxa 

8.11 

(±3.63) 

13.71 

3.70 

0.45 

% scraper  taxa 

14.41 
(±  2.69) 

7.55 

2.75 

0.19 

% eg  and  fg  taxa 

37.88 

(±2.25) 

5.27 

2.30 

0.06 

% predator  taxa 

18.09 

(±5.25) 

28.73 

5.36 

0.30 

% dominant  taxa 

15.71 

(±3.55) 

13.11 

3.62 

0.23 

HBP 

1.63 

(±0.26) 

0.07 

0.27 

0.17 

% asyncronous* 

0.06 

(±_0.02) 

0.00 

0.02 

0.30 

% semi-voitine* 

0.06 

(±0.02) 

0.00 

0.02 

0.26 

% uni-voltine^ 

0.57 

(±0.04) 

0.00 

0.04 

0.06 

Shannon  Diversity 
Index 

2.80 

(±) 

0.01 

0.05 

0.18 

% EPT  richness 

71.26 

(±5.69) 

0.34 

5.81 

0.08 

'Taxa  richness  is  the  number  of  taxa  per  sample;  ^EPT  richness  is  the  number  of  Ephemeroptera, 
Plecoptera,  and  Trichoptera  taxa  in  a sample;  ^cg  and  fg  is  the  number  of  collector-gatherer  and  filter- 
gatherer  taxa;^  HBI  is  Hilsenhoff  Biotic  Index;  ^ % voltinism  metrics  were  determined  to  be  too  variable 
for  further  consideration. 
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Table  13.  Mean  (4i  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Bozeman  Creek,  October  21,  1994  (N  = 3). 


metric 

mean 

(±  95%  C.I) 

variance 

s.  d. 

CV 

# of  organisms 

649.00 
C±  387.43) 

11721.10 

342.36 

0.53 

taxa  richness* 

20.67 
(+ ) 

8.33 

2.89 

0.14 

EPT  richness^ 

15.67 
C±  2.35) 

4.33 

2.08 

0.13 

# of  Plecoptera 
taxa 

4.00 

(±0.00) 

0.00 

0.00 

0.00 

#of 

Ephmeroptera 

taxa 

6.00 

(±1.13) 

1.00 

1.00 

0.17 

it  of  Trichoptera 
taxa 

5.67 

(±  1.30) 

1.33 

1.15 

0.20 

it  of  shredder  taxa 

4.00 

(±1.13) 

1.00 

1.00 

0.25 

it  of  scraper  taxa 

4.00 

(±1.13) 

1.00 

1.00 

0.25 

# of  eg  + fg  taxa-* 

4.33 

(±1.79) 

2.33 

1.58 

0.36 

it  of  predator  taxa 

4.00 

(±1.13) 

1.00 

1.00 

0.25 

% shredder  taxa 

19.22 
(±  3.37) 

8.85 

2.98 

0.16 

% scraper 
taxa 

19.22 

(±3.37) 

8.85 

2.98 

0.16 

% eg  + fg  taxa 

20.61 

(±5.23) 

21.35 

4.62 

0.22 

% predator  taxa 

19.22 

(±3.37) 

8.85 

2.98 

0.16 

% dominant  taxa 

20.75 

(±3.04) 

7.22 

2.69 

0.13 

iffiP 

1.31 

(±0.29) 

0.11 

0.34 

0.26 

% asyncronous* 

0.01 

(±0.03) 

0.00 

0.02 

1.73 

% semi-voltine* 

0.01 

(±  0.03) 

0.00 

0.02 

1.73 

% univoltine^ 

0.067 

(±0.14) 

0.02 

0.12 

0.19 

Shannon  Diversity 
Index 

2.31 

(±0.2) 

0.01 

0.86 

0.37 

%EPT  Richness 

75.88 

(±3.10) 

0.08 

2.74 

0.04 

'Taxa  richness:  number  of  taxa  per  sample;  ‘EPT  richness:  number  of  Ephemeroptera, 

^lecoptera,  and 

Trichoptera  taxa;  ^cg  and  fg:  number  of  collector-gatherer  and  filter-gatherer  taxa;  * HBI:Hilsenhoff 
Biotic  Index;^  % voltinism  metrics  were  determined  to  be  to  variable  for  further  consideration 
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Table  14.  Mean  (+.  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Bozeman  Creek,  February  5,  1995  (N  = 3). 


metric 

mean 

r-f  95^7.  C.I.) 

variance 

s.  d. 

CV 

# of  organisms 

259.67 

C±44.48) 

1545.30 

39.31 

0.15 

taxa  richness' 

26.33 

(±1.73) 

2.33 

1.53 

0.06 

EPT  richness^ 

19.CX) 

(±1.13) 

1.00 

1.00 

0.05 

it  of  Plecoptera 
taxa 

6.33 

(±1.73) 

2.33 

1.53 

0.24 

#of 

Ephmeroptera 

taxa 

6.67 

(±1.73) 

2.33 

1.53 

0.23 

it  of  Trichoptera 
taxa 

6.00 

(±0.00) 

0.00 

0.00 

0.00 

tt  of  shredder  taxa 

4.67 

(±  0.66) 

0.33 

0.58 

0.12 

it  of  scraper  taxa 

4.00 

(±1.13) 

1.00 

1.00 

0.25 

it  of  cg+fg  taxa^ 

5.33 

(±0.66) 

0.33 

0.58 

0.11 

it  of  predator  taxa 

6.67 

(±1.73) 

2.33 

1.53 

0.23 

% shredder  taxa 

17.84 

(±3.51) 

9.60 

3.10 

0.17 

% scraper 
taxa 

15.31 

(±4.87) 

18.50 

4.30 

0.28 

% eg  + fg  taxa^ 

20.22 

(±1.26) 

1.23 

1.11 

0.05 

% predator  taxa 

25.26 

(±6.10) 

29.05 

5.39 

0.21 

% dominant  taxa 

19.45 

(±5.98) 

27.93 

5.28 

0.27 

IffiT 

0.72 

(±0.14) 

0.33 

0.12 

0.17 

% asyncronous^ 

0.00 

(±0) 

0.00 

0.00 

0.00 

% semi-voltine^ 

0.05 

(±0.02) 

0.00 

0.02 

0.37 

% univoltine^ 

0.54 

(±0.033) 

0.00 

0.03 

0.05 

Shannon  Diversity’ 
Index 

2.69 

(±  0.20) 

0.15 

0.12 

0.05 

’Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 


* HBI  is  Hilsenhoff  Biotic  Index 

^ % voltinism  metrics  were  determined  to  be  too  variable  for  further  consideration 
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Table  15.  Mean  C±  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Bozeman  Creek,  May  5,  1995  (N  = 3). 


metric 

mean 

(+  95%  C.I.) 

variance 

s.  d. 

CV 

it  of  organisms 

396.67 
C±  102.30) 

8176.33 

90.40 

0.23 

taxa  richness* 

32.67 
(±  3.97) 

12.33 

3.51 

0.11 

EPT  richness^ 

25.00 

(±3.92) 

12.00 

3.46 

0.14 

it  of  Plecoptera 
taxa 

6.67 

(±  0.66) 

0.33 

0.58 

0.09 

#of 

Ephmeroptera 

taxa 

9.67 

(±2.61) 

2.33 

2.00 

0.11 

it  of  Trichoptera 
taxa 

8.33 

(±0.66) 

4.33 

0.58 

0.07 

a of  shredder  taxa 

3.33 

(±  0.66) 

0.33 

0.58 

0.17 

it  of  scraper  taxa 

3.33 

(±1.30) 

1.33 

1.15 

0.35 

a of  eg  and  fg 
taxa^ 

8.33 

(+1.73) 

2.33 

1.53 

0.18 

it  of  predator  taxa 

8.67 

(±  3.27) 

8.33 

2.89 

0.33 

% shredder  taxa 

10.27 

(±2.15) 

3.60 

1.90 

0.19 

% scraper 
taxa 

10.04 

(±3.13) 

7.67 

2.77 

0.28 

% eg  + fg  taxa^ 

26.05 
(±  8.69) 

59.01 

7.68 

0.29 

% predator  taxa 

26.23 

(±7.16) 

39.99 

6.32 

0.24 

% dominant  taxa 

6.67 

(±1.73) 

2.33 

1.53 

0.23 

IfflP 

1.52 

(±0.63) 

0.31 

0.56 

0.37 

% asyncronous^ 

0.05 

(±0.02) 

0.00 

0.01 

0.28 

% semi-voltine^ 

0.06 

(±0.06) 

0.00 

0.05 

0.85 

% univoltine^ 

0.50 

(±0.05) 

0.00 

0.04 

0.09 

Shannon  Diversity 
Index 

2.78 

(±0.3) 

0.19 

0.14 

0.05 

’Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 

^ % voltinism  metrics  were  determined  to  be  too  variable  for  further  consideration 
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Table  16.  Mean  (+.95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Bozeman  Creek,  September  20,  1995  (N=3). 


metric 

mean 

var. 

s.  d. 

CV 

# of  organisms 

673 

C±  257.45) 

51753.00 

227.50 

0.34 

taxa  richness' 

32.00 

(+1.13) 

1.00 

1.00 

0.03 

EPT  richness^ 

24.33 

(+0.66) 

0.33 

0.58 

0.02 

M of  Plecoptera 
taxa 

5.33 

(+0.66) 

0.33 

0.58 

0.11 

#of 

Ephemeroptera 

taxa 

9.33 

C+  0.66) 

0.33 

0.58 

0.06 

it  of  Trichoptera 
taxa 

9.67 

(±1.73) 

2.33 

1.53 

0.16 

tt  of  shredder 
taxa 

4.33 

(±1.30) 

1.33 

1.15 

0.27 

tt  of  scraper  taxa 

4.00 

(±1.00) 

1.00 

1.00 

0.25 

tt  of  cg+fg  taxa^ 

5.33 

(±  0.66) 

0.33 

0.58 

0.11 

tt  of  predator 
taxa 

6.67 

(±  1.73) 

2.33 

1.53 

0.23 

% shredder 
taxa 

13.49 

(±3.74)" 

10.93 

3.31 

0.25 

% scraper  taxa 

12.48 

(±3.29) 

8.46 

2.91 

0.23 

% cg-»-fg  taxa^ 

16.68 

(±2.11) 

3.46 

1.86 

0.11 

% predator  taxa 

20.75 

(±4.72) 

17.4 

4.17 

0.20 

% dominant  taxa 

12.67 
(±  1.05 

0.87 

0.93 

0.07 

ifflp 

1.61 

(±  0.24) 

0.04 

0.21 

0.13 

Shannon 
diversity  index 

3.00 
(±  0.30) 

0.00 

0.07 

0.00 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Two  richness  indices,  four  diversity  indices,  and  five  evenness  indices  were 
computed  for  the  five  sample  times  from  Bozeman  Creek.  Results  are  presented  in 
Table  16.  Shannon  diversity  index  is  the  most  commonly  used  diversity  index  by  1) 
ecologists  (Ludwig  and  Reynolds  1988),  2)  most  rapid  assessment  programs  (Barbour 
et  al.  1995),  and  3)  is  recommended  for  use  by  the  Water  Quality  Division  of  Montana 
(Bukantis  1995b).  Further  analysis  of  Shannon  diversity  index  was  made  for  seasonal 
variability  in  Bozeman  Creek. 

The  mean,  variance,  s.d.  and  CV  of  the  Shannon  diversity  index  values  are 
presented  in  Table  18.  Results  of  this  analysis  suggested  that  there  was  little  variability 
in  species  diversity  between  each  riffle  sample  within  each  season  and  between  each 
season  as  shown  by  the  small  CV  values. 

To  a limited  extent,  the  percent  dominant  taxa  metric  can  also  be  used  to 
analyze  species  diversity  (evenness).  The  % dominant  taxa  as  measured  by  relative 
abundance,  comprised  only  6.67%  C±.  1.73%)  of  the  organisms  in  the  Spring  1995 
samples  (Appendix  21).  On  the  other  hand,  20.61%  (+.  3.04%)  and  19.45%  (+. 
5.98%)  of  the  organisms  in  the  Autumn  1994  and  Winter  1995  samples  were  from  one 
taxon,  the  Plecopteran  shredder  Zapada  spp.  Summer  1994  and  Summer  1995,  % 
dominant  taxon  accounted  for  15.71%  (i  3.55%)  and  12.67%  (+.  1.05%)  of  the  total 


number  of  organisms,  respectively. 
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Table  17.  Species  diversity  indices  for  sixteen  samples  from  five  sampling  periods  in  Bozeman  Creek,  Spring  1995, 
Summer  1994  (Sml,  Smr2,  Smr3,and  Smr4)  , Fall  1994,  Winter  1995,  and  Summer  1995  (Smrl,Smr2,  and  Smr3). 
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Table  18.  Mean,  variance,  standard  deviation  (s.d.),  and  coefficient  of  variation 
(CV)  for  Shannon  diversity  index  for  four  seasons  (including  two  summer 
seasons)  from  Bozeman  Creek,  1994  and  1995. 


Summer  1994 
(N=4) 

Fall  1994 
(N=3) 

Winter  1994 
(N=3) 

Spring  1995 
(N=3) 

Summer 
1995 
(N  = 3) 

mean 

2.80 

2.31 

2.69 

2.78 

3.00 

variance 

0.00 

0.01 

0.02 

0.02 

0.00 

s.d. 

0.05 

0.09 

0.12 

0.14 

0.07 

CV 

0.02 

0.04 

0.05 

0.05 

0.0 

Samples  collected  in  late  summer  1995  from  Bozeman  Creek  probably  reflect 
the  transition  between  early  summer  macroinvertebrate  community  composition  ( 
summer  1994  samples)  and  between  the  autumn  macroinvertebrate  community 
composition  (autumn  1994  samples).  The  number  of  organisms  collected  in  late 
summer  1995  was  673.00  (+.  257.45)  which  was  closer  to  the  autumn  1994  number  of 
organisms  mean  of  649.00  C±  387.43),  while  the  mean  number  of  organisms  collected 
in  early  summer  1994  was  315.75  (±_  63.56).  Number  of  taxa  (taxa  richness)  for  late 
summer  1995  was  32.0  (+.1.13)  which  is  more  similar  to  the  mean  number  of  taxa  for 
the  early  summer  1994  collections  at  27.75C+  1.24).  The  autumn  1994  mean  number 
of  taxa  was  only  20.67  (+.  3.27). 

A quantitative  similarity  index  (QSI)  for  taxa  was  used  to  quantify  the 
similarities  between  samples  within  each  season  and  between  seasons  (Table  19).  The 
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QSI  used  % relative  abundances  of  each  species  from  each  sample  and  compared  this 
percentage  between  all  three  samples  to  give  a total  % similarity  between  samples  for 
each  season.  For  Bozeman  Creek,  all  three  samples  were  used  from  each  season  to 
determine  a QSI  value,  except  for  Summer  1994.  Since  Summer  1994  had  four 
samples,  a mean  value  of  the  four  combinations  of  three  samples  collected  in  Summer 
1994  was  used  (Table  20). 


Table  19.  Quantitative  Similarity  Index  values  between  the  three  samples  for  each 
season.  Summer  1994,  Autumn  1994,  Winter  1995,  Spring  1995,  and  Summer  1995 
including  mean  (+.  95%  C.  I.),  variance  (var.),  s.  d.  and  CV. 


smmr 

94' 

atmn 

94 

wntr 

95 

spmg 

95 

smmr 

95 

mean 

var. 

s.d. 

CV 

QSI 

43.05 

45.37 

52.48 

58.82 

69.92 

53.93 

±9.53 

10.76 

10.87 

0.3 

‘Summer  1994  QSI  was  from  the  mean  value  from 


bur  combinations  of  three  samples 


used  to  standardize  each  season  at  three  samples  (Table  20). 


Table  20.  Quantitative  Similarity  Index  (QSI)  values  for  four  combinations  of  three 
samples  (1,2, 3,4)  from  Bozeman  Creek  summer  1994  including  mean  (+.95%  C.I.), 
variance,  s.  d.,  and  CV. 


1,2,3 

1,3,4 

2,3,4 

1,2,4 

mean 

variance 

s.d. 

CV 

QSI 

48.32 

41.19 

44.69 

37.99 

43.05 

(± 

4.36) 

19.84 

4.45 

0.10 

QSI  values  used  for  Bozeman  Creek  seasonal  variability  showed  that  Summer 
1994  and  Spring  1995  samples  were  least  similar  (Table  21).  The  QSI  value  was  0.334 
or  33.4%  similarity.  The  most  similar  seasons  using  QSI  were  Winter  1995  and 
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Spring  1995  with  a value  of  0.82  (82%  similar).  The  only  other  seasons  greater  than 
50%  similarity  were  Fall  1994  and  Winter  1995  with  a QSI  value  of  0.516  (Table  21). 


Table  21.  Quantitative  Similarity  Index  (QSI)  comparisons  between  mean  % relative 
abundances  for  Summer  1994,  Autumn  1994,  Winter  1995,  and  Spring  1995  from 
Bozeman  Creek  macroinvertebrate  samples.  


Season 

QSI 

Summer  94/Spring  95 

33.40 

Autumn  94/Spring  95 

34.50 

Summer  94/Autumn  94 

35.13 

Winter  95/Summer  94 

38.67 

Autumn  94/Winter  95 

51.61 

Winter  95/Spring  95 

82.00 

The  modified  Hilsenhoff  Biotic  Index  (HBI)  is  a measure  of  a taxons' 
sensitivity  to  organic  pollution.  HBI  values  range  from  0 for  least  tolerant  taxa  to  10 
for  most  tolerant  taxa.  HBI  values  were  lowest  for  Winter  1995  samples  with  a mean 
of  0.72,  indicating  that  the  most  sensitive  taxa  to  organic  pollution  occurred  in  the 
benthos  during  winter  season. 

Mean  and  95%  Cl's  for  seventeen  metrics  for  the  five  seasons  sampled  in 
Bozeman  Creek  are  in  Figures  23  through  39. 
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Figure  23.  Mean  {±95%  Cl)  number  of  organisms  for  Bozeman  Creek. 
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Figure  24.  Mean  C±  95%  Cl)  number  of  taxa  (taxa  richness)  for  Bozeman  Creek. 
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Figure  25.  Mean  C±  95%  Cl)  number  of  EPT  taxa  for  Bozeman  Creek. 
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Figure  26.  Mean  (±95%  Cl)  number  of  Ephemeroptera  taxa  for  Bozeman  Creek. 
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Figure  27.  Mean  (±95%  Cl)  number  of  Plecoptera  taxa  for  Bozeman  Creek. 
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Figure  28.  Mean  C±  95%  Cl)  number  of  Trichoptera  taxa  for  Bozeman  Creek. 
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Figure  29.  Mean  (±95%  Cl)  number  of  shredder  taxa  for  Bozeman  Creek. 
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Figure  30.  Mean  (±95%  Cl)  number  of  scraper  taxa  for  Bozeman  Creek. 
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Figure  31.  Mean  (±95%  Cl)  number  of  collector-gatherer  and  filter-gatherer  taxa  for  Bozeman  Creek. 
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Figure  32.  Mean  C±  95%  Cl)  number  of  predator  taxa  for  Bozeman  Creek. 
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Figure  33.  Mean  (Ji  95%  Cl)  % shredder  taxa  for  Bozeman  Creek. 
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Figure  34.  Mean  C±.  95%  Cl)  % scraper  taxa  for  Bozeman  Creek. 
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Figure  35.  Mean  C±  95%  Cl)  % collector-gatherer  and  Filter-gatherer  taxa  for  Bozeman  Creek. 
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Figure  36.  Mean  (i  95%  Cl)  % predator  taxa  for  Bozeman  Creek. 
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Figure  37.  Mean  (±95%  Cl)  % dominant  laxa  for  Bozeman  Creek. 
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Figure  38.  Mean  (±95%  Cl)  HBI  values  for  Bozeman  Creek. 
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Figure  39.  Mean  {±,95%  Cl)  Shannon  Diversity  Index  values  for  Bozeman  Creek. 
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The  Three  Bear  Creeks  1994  and  1995 

Principal  Components  Analysis  comparing  raw  abundance  data  from  each 
sample  as  a sampling  unit  for  the  three  Bear  Creeks  (N.  F.  Bear  Creek,  M.  F.  Bear 
Creek,  and  Bear  Creek)  suggested  that  Bear  Creek,  1994  (N  = l)  was  less  similar  to  any 
of  the  other  samples  collected  from  all  of  the  three  Bear  Creeks  (Table  22  and  Figure 
40).  Bear  Creek,  1995  samples  (N  = 3)  formed  a looser  cluster  than  did  N.  F.  Bear 
Creek  1995  and  M.  F.  Bear  Creek  1994  and  1995  (Table  22  and  Figure  40).  N.  F. 
Bear  Creek  1994  (N  = 1)  also  clustered  separately. 


Table  22.  Eigenvectors  defining  the  first  three  principal  components  between  sampling 
units  for  N.  F.  Bear  Cr.,  M.  F.  Bear  Cr.  and  Bear  Creek,  1994  and  1995. 


Sampling  Unit 

I 

n 

m 

N.F.BearCr.1995 

1 

-0.745 

0.132 

-0.983 

2 

-0.951 

0.479 

-1.374 

3 

-0.334 

-0.710 

-0.268 

N.F.Bear  Cr.l994 

4 

-0.085 

1.220 

0.162 

M.F.  Bear  Cr.  1995 

5 

-0.259 

0.123 

0.369 

6 

-0.520 

-0.097 

0.355 

7 

-0.458 

-0.518 

1.181 

M.F.  Bear  Cr.  1994 

8 

-1.009 

-0.015 

1.124 

Bear  Cr.  1995 

9 

0.008 

-0.346 

-0.421 

10 

1.724 

-1.454 

-0.303 

11 

0.806 

-0.488 

-0.078 

Bear  Cr.  1994 

12 

1.824 

1.674 

0.235 

Mean,  95%  Cl,  variance,  s.d.,  and  CV  for  metrics  calculated  for  the  three  Bear 
Creeks  are  in  Tables  23  through  26. 
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Figure  40.  Principal  Componenls  Analysis  for  the  three  Bear  Creeks. 
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Table  23.  Mean  (i  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  N.  F.  Bear  Creek,  1995  (N  = 3). 


metric 

mean 

(+  95%  C.  I.) 

variance 

s.  d. 

CV 

# of  organisms 

375.33 
C±  120.30) 

11301.33 

106.31 

0.28 

taxa  richness' 

23.0 

(+3.39) 

9.00 

3.00 

0.13 

EPT  richness^ 

19.33 
(+  2.36) 

4.33 

2.08 

0.11 

# of  Plecoptera 
taxa 

3.00 

(+1.13) 

1.00 

1.00 

0.33 

#of 

Ephemeroptera 

taxa 

10.00 

(+0.00) 

0.00 

0.00 

0.00 

# of  Trichoptera 
taxa 

6.33 

(+1.73) 

2.33 

1.53 

0.24 

# of  shredder 
taxa 

2.00 

(+0.00) 

0.00 

0.00 

0.00 

M of  scraper  taxa 

2.00 

C±0.00) 

0.00 

0.00 

0.00 

# of  cg+fg  taxa"* 

8.67 

(+0.65) 

0.33 

0.58 

0.07 

H of  predator 
taxa 

5.67 

(+  1.73) 

2.33 

1.53 

0.27 

% shredder 
taxa 

8.80 

(+1.31) 

0.01 

1.16 

0.13 

% scraper  taxa 

8.80 

(+1.31) 

0.01 

1.16 

0.13 

% cg+fg  taxa^ 

38.30 
(+  8.09) 

0.51 

7.15 

0.19 

% predator  taxa 

24.50 
(+  6.07) 

0.29 

5.36 

0.22 

% dominant  taxa 

50.30 
(±  20.30) 

3.20 

17.90 

0.36 

HBP 

1.87 

(±  0.46) 

0.16 

0.40 

0.21 

Shannon 
diversity  index 

1.91 

(+0.54) 

0.23 

0.48 

0.25 

%EPT  by  taxa 

84.41 

(±8.16) 

0.51 

7.21 

0.09 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  24.  Mean  (+_  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  M.  F.  Bear  Creek,  1995  (N  = 3). 


metric 

mean 

variance 

s.  d. 

CV 

M of  organisms 

361.33 
C±  71.00) 

3936.33 

62.74 

0.17 

taxa  richness* 

20.00 

(+4.53) 

16.00 

4.00 

0.2 

EPT  richness^ 

17.33 

(+2.36) 

4.33 

2.08 

0.12 

a of  Plecoptera 
taxa 

2.67 

C+ 0.65) 

0.33 

0.58 

0.22 

#of 

Ephemeroptera 

taxa 

8.33 

(+  0.65) 

0.33 

0.58 

0.07 

M of  Trichoptera 
taxa 

6.33 

(+1.31) 

1.33 

1.15 

0.18 

it  of  shredder 
taxa 

2.00 

(±0.00) 

0.00 

0.00 

0.00 

it  of  scraper  taxa 

2.67 

(±0.65) 

0.33 

0.58 

0.22 

a of  cg+fg  taxa-* 

7.33 

(±0.65) 

0.33 

0.58 

0.08 

it  of  predator 
taxa 

4.67 

(±0.65) 

0.33 

0.58 

0.12 

% shredder 
taxa 

10.28 
(±  2.37) 

0.04 

2.1 

0.02 

% scraper  taxa 

13.33 

(±1.63) 

0.02 

1.44 

0.11 

% cg+fg  taxa^ 

37.36 

(±6.33) 

0.31 

5.59 

0.15 

% predator  taxa 

24.03 
(±  7.09) 

0.39 

6.27 

0.26 

% dominant  taxa 

51.00 

(±5.67) 

0.25 

5.02 

0.1 

ifflP 

1.81 

(±0.08) 

0.01 

0.07 

0.04 

Shannon 
diversity  index 

1.85 

(±  0.29) 

0.06 

0.25 

0.13 

%EPT  by  taxa 

87.64 
(±  8.57) 

0.57 

( 

7.57 

0.09 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  25.  Mean  C±  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Bear  Creek,  1995  (N  = 3). 


metric 

mean 

var. 

s.d. 

CV 

it  of  organisms 

286.00 
(±  126.31) 

12459.00 

111.62 

0.39 

taxa  richness* 

23.67 

(+1.73) 

2.33 

1.53 

0.06 

EPT  richness^ 

21.00 
(+  2.99) 

7.00 

2.65 

0.13 

it  of  Plecoptera 
taxa 

3.33 

(+1.73) 

2.33 

1.53 

0.46 

#of 

Ephemeroptera 

taxa 

10.67 

(+0.65) 

0.33 

0.58 

0.05 

it  of  Trichoptera 
taxa 

7.00 

(±1.13) 

1.00 

1.00 

0.14 

it  of  shredder 
taxa 

1.33 

(+0.65) 

0.33 

0.58 

0.43 

it  of  scraper  taxa 

3.00 

(+0.00) 

0.00 

0.00 

0.00 

a of  cg+fg  taxa-* 

10.00 

(±1.13) 

1.00 

1.00 

0.1 

it  of  predator 
taxa 

5.67 

(±  2.85) 

6.33 

2.52 

0.44 

% shredder 
taxa 

5.75 

(±3.27) 

0.08 

2.89 

0.5 

% scraper  taxa 

12.71 

(±0.95 

0.00 

0.84 

0.07 

% cg+fg  taxa-* 

42.19 

(±1.82) 

0.02 

1.61 

0.04 

% predator  taxa 

23.55 
(±  10.49) 

0.86 

9.27 

0.39 

% dominant  taxa 

20.19 

(±3.49) 

0.01 

3.08 

0.15 

hbP 

1.02 

(±0.27) 

0.06 

0.24 

0.21 

Shannon 
diversity  index 

2.59 

(±0.10) 

0.01 

0.09 

0.03 

%EPT  by  taxa 

88.49 
(±  6.55) 

0.33 

5.78 

0.07 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 


S ..■*^7 


i'Ji .,  ' i’-  '■  iM 


ST. 


10  .nsti  \ “.? « ,i  i\  7 / 


-■  IS' ■' 


isl 


J 


■ 'J'j 


■ -'d 


i’ 


J,' 


•■  ij  . 

r i 1 W 


. ' ■' ■A  '-  ' 'P\  It-V" 


hw;^y 


WMMkx.^Jh- 


0 


1 


li  ' '■'  • - 1 


r>|»r.4i^.  r(t« 

l’^'  ...',,  . , % 


. • • »>ww.  «bv 

y .^j',  ► 


J if'*’ 

m-i' 


M 


i.Q 


1 #. 


1. 


- /»*♦«* 

J V^' 


■U,:.  - 


"L. 

1 _>*HW»ivi.» 


Mr 


'’Hil' 


^ ■ .^u— ,v 

f "M'  •. 


'-‘ri'"-  (?$' 

. ■: 

f ‘i^^i-*vjr»  ^ '.r'-ftffAfe-'  .!.j 

'Mm.  . ^ 


-iWr<  * -y^ 


,-tw 

■< " t 

r, 


' '?  ■ , . ■ . j 


- f V'rtf—  *ir— ■J.  *JI^‘  ’,'1^ 

» J:' 


>A/^  ■ 


'■V  ‘S'-;  ■ 


J[i7] 


j ^):0'  ^ ' ''  , ; 


5’  r - 


s^'; 


W: 


Aa)i'} 


- 


'in 


.a- 


w 


M, 


..y 


K.  M; 


■m 


■.;.%* 


■:vii. 


Li?« 


S; 


96 


Table  26.  Metric  values  for  macroinvertebrate  samples  collected  from  N.  F.  Bear 
Creek,  M.  F.  Bear  Creek,  and  Bear  Creek  in  1994.  One  sample  was  collected  from 
each  creek. 


metric 

N.F.  Bear 
Creek  1994 

M.F.  Bear 
Creek  1994 

Bear  Creek 
1994 

# of  organisms 

383.00 

233.00 

408.00 

taxa  richness' 

22.00 

20.00 

25.00 

EPT  richness^ 

17.00 

15.00 

19.00 

# of  Plecoptera  taxa 

3.00 

2.00 

4.00 

it  of  Ephemeroptera  taxa 

9.00 

7.00 

8.00 

tf  of  Trichoptera  taxa 

5.00 

6.00 

7.00 

tt  of  shredder  taxa 

2.00 

0.00 

1.00 

tt  of  scraper  taxa 

2.00 

2.00 

4.00 

ti  of  cg+fg  taxa^ 

8.00 

7.00 

6.00 

it  of  predator  taxa 

5.00 

4.00 

5.00 

% shredder  taxa 

9.09 

0.00 

4.00 

% scraper  taxa 

9.09 

10.00 

16.00 

% cg+fg  taxa-* 

36.36 

35.00 

24.00 

% predator  taxa 

22.73 

20.00 

2.00 

% dominant  taxa 

24.3 

41.44 

29.66 

hdbP 

1.14 

1.49 

0.91 

Shannon  diversity  index 

2.18 

1.96 

2.09 

%EPT^  by  taxa 

77.27 

75.00 

76.00 

'Taxa  richness  is  the  number  of  laxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 


Impaired  vs.  Unimpaired  (Reference)  Streams 

Box  and  whisker  plots  using  the  median,  minimum,  maximum  and  25th  to  75th 
percentile  values  comparing  the  seven  reference  streams  (N  = 34  samples)  and  three 
most  severely  impaired  streams.  Soda  Butte  Creek  (site  2),  Divide  Creek,  and  Hyalite 
Creek  (N  = 12  samples)  suggested  that  out  of  the  19  metrics  compared,  5 metrics  were 
probably  not  useful  for  detecting  impairment  based  on  25th  to  75th  percentile  overlap 
(Figure  41  through  Figure  59).  These  were  %EPT  taxa,  % dominant  taxa,  modified 
HBI,  % Chironomidae  organisms,  and  % collector-gatherer  and  filter-gatherer  taxa.  T- 
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test  comparisons  of  the  number  of  organisms  (t=4.12;  p=0.001;  df=9),  number  of 
Ephemeroptera  taxa  (t=3.96;  p=0.002;  df=9),  % dominant  taxa  (t=0.594;  p=0.285; 
df=9)),  and  % Chironimidae  organisms  (t=0.95;  p= 0.285;  df=9)  between  reference 
and  impaired  sites  using  the  mean  of  combined  samples  from  each  stream  (N=7 
reference;  N=3  impaired)  showed  that  the  box  and  whisker  plot  method  was  reliable. 

Divide  Creek 

Mean,  95%  Cl,  variance,  s.d.,  and  CV  for  Divide  Creek  1995  samples  (N  = 3) 
are  in  Table  27.  Metric  scores  from  the  Divide  Creek,  1994  sample  is  in  Table  28. 

Soda  Butte  Creek 

Principal  components  analysis  (PCA)  of  macroinvertebrates  from  Soda  Butte 
Creek,  1995  showed  that  macroinvertebrate  communities  grouped  into  three  distinct 
faunal  units.  These  were  1)  the  macroinvertebrate  community  above  the  McLaren 
tailing  complex,  2)  sites  2 Km  and  4 Km  directly  below  the  McLaren  tailing  complex 
and  3)  sites  8 Km,  11  Km,  and  16  Km  which  were  further  downstream  from  the 
tailings  and  were  less  influenced  by  the  tailings  due  to  dilution  affects  (Figure  60  and 
Table  29).  PCA  analysis  of  Pebble  Creek,  which  is  adjacent  to  Soda  Butte  Creek  and 
considered  'unimpaired'  grouped  macroinvertebrates  into  only  two  faunal  units  (Figure 


2). 
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Figure  41.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  number  of  organisms. 
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ref  imp 

min..max.,  25-75th  percentile,  median 


Figure  42.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  number  of  taxa. 
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Figure  43.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25th  to  75th  percentile  values  for  number  of  EPT  taxa. 
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Figure  44.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75di  percentile  values  for  number  of  Ephemeroptera  taxa. 
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Figure  45.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25 Ui  to  75tli  percentile  values  for  number  of  Plecoptera  taxa. 
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Figure  46.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25Ui  to  75tli  percentile  values  for  number  of  Trichoptera  taxa. 
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Figure  47.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25th  to  75th  percentile  values  for  number  of  shredder  taxa. 


104 


min.,max.,25-75th  percentile,  median 


Figure  48.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  number  of  scraper  taxa. 
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Figure  49.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25th  to  75th  percentile  values  for  number  of  collector-gatherer  and  filter-gatherer  laxa. 
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Figure  50.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25Ui  to  75lh  percentile  values  for  number  of  predator. 
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Figure  51.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  % shredder  taxa. 
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Figure  52.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25th  to  75th  percentile  values  for  % scraper  taxa 
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Figure  53.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25th  to  75Ui  percentile  values  for  % collector-gatherer  and  filter-gatherer  laxa. 
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Figure  54.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  % predator  taxa. 
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Figure  55.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75th  percentile  values  for  % dominant  taxa. 
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Figure  57.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  % Chironomidae  organisms. 
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Figure  58.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25tli  to  75tli  percentile  values  for  HBI  values. 
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Figure  56.  Impaired  vs.  unimpaired  (reference)  stream  box  and  whisker  plots  showing  median,  minimum, 
maximum,  and  25th  to  75th  percentile  values  for  % EFT  taxa. 
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Table  27.  Mean  (±95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Divide  Creek,  1995  (N  = 3). 


metric 

mean 

(±  95%  C.I.) 

variance 

s.d. 

CV 

M of  organisms 

65.67 
(±  17.02) 

226.33 

15.04 

0.23 

taxa  richness* 

11.00 

(+1.13) 

1.00 

1.00 

0.09 

EPT  richness^ 

8.33 

(+1.73) 

2.33 

1.53 

0.18 

a of  Plecoptera 
taica 

2.00 

(+0.00) 

0.00 

0.00 

0.00 

#of 

Ephemeroptera 

taxa 

6.33 

(+1.73) 

2.33 

1.53 

0.24 

# of  Trichoptera 
taxa 

0.00 

C±0.00) 

0.00 

0.00 

0.00 

# of  shredder  taxa 

1.00 

(+0.00) 

0.00 

0.00 

0.00 

# of  scrapier  taxa 

0.67 

(+0.65) 

0.33 

0.58 

0.87 

M of  eg  and  fg 
taxa’ 

5.67 

(+1.31) 

1.33 

1.15 

0.2 

M of  predator  taxa 

1.00 

(+0.00) 

0.00 

0.00 

0.00 

% shredder  taxa 

9.14 

(±0.94) 

0.00 

0.83 

0.09 

% scraper  taxa 

6.11 

(+  6.06) 

0.29 

5.36 

0.87 

% eg  + fg  taxa’ 

51.26 
(±  7.39) 

0.43 

6.53 

0.13 

% predator  taxa 

9.14 

(+0.94) 

0.00 

0.83 

0.09 

% dominant  taxa 

40.27 
(+  6.36) 

0.32 

5.62 

0.14 

ifflp 

0.78 

(+0.31) 

0.07 

0.27 

0.35 

Shannon  Diversity 
Index 

1.84 

(+0.12) 

0.01 

0.11 

0.06 

%EPT  richness^ 

75.66 
(±  11.93) 

1.11 

10.54 

0.14 

‘Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  28.  Metric  scores  from  Divide  Creek,  July  16,  1994  (one  sample). 


metric 

number 

# of  organisms 

105 

taxa  richness* 

9 

EPT  richness^ 

8 

# of  Plecoptera  taxa 

1 

# of  Ephemeroptera 

7 

# of  Trichoptera  taxa 

0 

# of  shredder  taxa 

0 

# of  scraper  taxa 

1 

# of  eg  and  fg  taxa* 

7 

# of  predator  taxa 

1 

% shredder  taxa 

0 

% scraper  taxa 

11.11 

%cg  and  fg  taxa^ 

77.77 

% predator  taxa 

11.11 

% dominant  taxa 

29.52 

HBP 

0.70 

Shannon  Diversity  Index 

1.76 

%EPT  richness^ 

88.89 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  29.  Eigenvectors  defining  the  first  three  principal  components  between  18 
sampling  units  from  Soda  Butte  Creek,  1995. 

Sample  Unit I H HI 


Km  1 

1 

2 

3 

-0.61 

-0.56 

-0.67 

1.30 

1.67 

0.69 

-0.39 

-0.27 

0.17 

Km2 

4 

-0.88 

-0.53 

0.04 

5 

-0.86 

-0.35 

0.13 

6 

-0.78 

-0.58 

0.04 

Km4 

7 

-0.62 

-0.52 

0.02 

8 

-0.56 

-0.59 

0.00 

9 

-0.36 

-0.71 

-0.04 

Km8 

10 

0.15 

-0.11 

0.75 

11 

0.14 

-0.13 

0.33 

12 

1.10 

0.69 

1.69 

Km  11 

13 

1.12 

-0.49 

-0.46 

14 

0.93 

-0.16 

0.02 

15 

0.80 

-0.27 

-0.31 

Km  16 

16 

0.41 

-0.36 

-0.28 

17 

0.32 

0.03 

-0.12 

18 

0.94 

0.41 

-1.31 

Mean,  95%  Cl,  variance,  s.d.,  and  CV  for  metrics  calculated  for  Soda  Butte 
Creek  km  1 and  km  2,  1995,  are  in  Table  30  and  Table  31.  Metric  values  calculated 
for  Soda  Butte  Creek  km  1 and  km  2,  1994  are  in  Table  32.  Mean  and  95%  Cl 
longitudinal  and  impairment  metric  comparisons  are  in  Figures  61  through  76. 

Reference  streams.  Cottonwood  Creek.  Hyalite  Creek,  and  Squaw  Creek 

Mean,  95%  Cl,  variance,  s.d.,  and  CV  for  the  seven  reference  streams 
combined  are  in  Table  33.  Mean,  95%  Cl,  variance,  s.d.,  and  CV  for  Cottonwood 
Creek,  Hyalite  Creek,  and  Squaw  Creek  are  in  Tables  34  through  36. 
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Figure  60.  Principal  Components  Analysis  for  Soda  Butte  Creek,  1995 
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Table  30.  Mean  (±_  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Soda  Butte  Creek  site  km  1,  1995  (N  = 3). 


metric 

mean 

(±  95%  C.I.) 

variance 

s.d. 

CV 

M of  organisms 

238.33 

(±62.21) 

3022.33 

54.98 

0.23 

taxa  richness^ 

22.67 
(+  2.36) 

4.33 

2.08 

0.09 

EPT  richness^ 

17.00 

(±1.96) 

3.00 

1.73 

0.10 

M of  Plecoptera 
taxa 

3.67 

(±1.73) 

2.33 

1.53 

0.42 

Ephemeroptera 

taxa 

6.67 

(±0.65) 

0.33 

0.58 

0.09 

# of  Trichoptera 
taxa 

6.67 

(±0.65) 

0.33 

0.58 

0.09 

it  of  shredder  taxa 

2.00 

(±0.00) 

0.00 

0.00 

0.00 

it  of  scraper  taxa 

1.00 

(±0.00) 

0.00 

0.00 

0.00 

# of  eg  and  fg 
taxa^ 

5.67 

(±0.65) 

0.33 

0.58 

0.10 

it  of  predator  taxa 

7.00 

(±1.96) 

3.00 

1.73 

0.25 

% shredder  taxa 

9.69 

(±1.00) 

0.00 

0.50 

0.05 

% scraper  taxa 

4.85 

(±0.00) 

0.00 

0.30 

0.05 

% eg  + fg  taxa-* 

27.42 

(±3.00) 

0.06 

2.50 

0.09 

% predator  taxa 

33.79 

(±9.00) 

0.60 

7.82 

0.23 

% dominant  taxa 

24.72 

(±8.00) 

0.48 

6.98 

0.28 

ffflP 

1.91 

(±0.37) 

0.11 

0.33 

0.17 

Shannon  Diversity 
Index 

2.35 

(±  0.09) 

0.01 

0.08 

0.03 

%EPT  richness^ 

74.97 
(±  2.00) 

0.00 

1.90 

0.03 

‘Taxa  richness  is  the  number  of  taxa  p>er  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample. , 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  31.  Mean  (±_  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  Soda  Butte  Creek  site  km  2,  1995  (n=3). 


metric 

mean 

(+  95%  C.I.) 

variance 

s.d. 

CV 

# of  organisms 

19.67 

(+6.23) 

30.33 

5.51 

0.28 

taxa  richness* 

8.00 

(+1.96) 

3.00 

1.73 

0.22 

EPT  richness^ 

6.00 

(+1.96) 

3.00 

1.73 

0.29 

it  of  Plecoptera 
taxa 

1.00 

(+1.13) 

1.00 

1.00 

1.00 

#of 

Ephemeroptera 

taxa 

2.00 

(±0.00) 

0.00 

0.00 

0.00 

it  of  Trichoptera 
taxa 

3.00 

(+  2.26) 

4.00 

2.00 

0.67 

it  of  shredder  taxa 

0.67 

(±  0.65) 

0.33 

0.58 

0.87 

it  of  scraper  taxa 

0.00 

(+0.00) 

0.00 

0.00 

0.00 

it  of  eg  andfg 
taxa^ 

2.67 

(+0.65) 

0.33 

0.58 

0.22 

it  of  predator  taxa 

. 1.33 

(+1.73) 

2.33 

1.53 

1.15 

% shredder  taxa 

8.10 

(±  8.30) 

0.54 

7.30 

0.91 

% scraper  taxa 

0.00 

(±0.00) 

0.00 

0.00 

0.00 

% eg  + fg  taxa^ 

33.81 

(±8.90) 

0.60 

7.87 

0.23 

% predator  taxa 

14.76 
(±  16.98) 

2.25 

15.01 

1.01 

% dominant  taxa 

41.93 
(±  8.70) 

0.60 

7.70 

0.18 

imp 

1.66 

(±0.79) 

0.49 

0.70 

0.42 

Shannon  Diversity 
Index 

1.58 

(±  0.29) 

0.07 

0.26 

0.16 

%EPT  organisms 

63.18 
(±  17.93) 

2.51 

15.84 

0.25 

% Chironimidae 

36.80 
(±  17.93) 

2.51 

15.84 

0.43 

‘Taxa  richness  is  the  number  of  laxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  laxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  laxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  32.  Metrics  calculated  for  Soda  Butte  Creek  sites  km  1 and  km  2,  1994. 


metric 

Soda  Butte 
Creek  Site  1 

Soda  Butte 
Creek  Site  2 

# of  organisms 

354.00 

16.00 

taxa  richness’ 

23.00 

7.00 

EPT  richness^ 

15.00 

6.00 

a of  Plecoptera  taxa 

4.00 

2.00 

# of  Ephemeroptera 
taxa 

8.00 

3.00 

M of  Trichoptera  taxa 

3.00 

1.00 

# of  shredder  taxa 

2.00 

1.00 

# of  scraper  taxa 

2.00 

0.00 

a of  eg  and  fg  taxa-* 

7.00 

4.00 

# of  predator  taxa 

7.00 

1.00 

% shredder  taxa 

8.70 

14.29 

% scraper  taxa 

8.70 

0.00 

% eg  -1-  fg  taxa-’ 

30.43 

57.14 

% predator  taxa 

30.43 

14.29 

% dominant  taxa 

23.73 

25.00 

imp 

1.68 

1.50 

Shannon  Diversity  Index 

2.45 

1.79 

%EPT  richness^ 

0.65 

0.86 

'Taxa  richness  is  the  number  of  laxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Figure  61.  Mean  (95%  Cl)  number  of  organisms  for  Soda  Butte  Creek,  1995. 
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Figure  62.  Mean  (95%  Cl)  number  of  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  63.  Mean  (95%  Cl)  number  of  EPT  taxa  for  Soda  Butte  Creek,  1995. 


125 


10  %96  pue  ueeuj  1 


exB]  id3  jo  # 


Figure  64.  Mean  (95%  Cl)  number  of  Ephemeroptera  taxa  for  Soda  Butte  Creek,  1995. 


126 


10  %96  PUB  UB0LU 

: ; t 1 : 

i i i i : 

j : : I I 

■-4-1 — i 

t I i i : 


i-Mi — i 

I j 

i I 

i I 

! i 


I 

i 


O)  00  t^  CD  m Tt  CO 

BXBl  BJ0ldOJ0LU0Lld3  p # 


00 

E 

^ </> 
X CD 


E 


CN 


CO 


Figure  65.  Mean  (95%  Cl)  number  of  Plecoptera  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  66.  Mean  (95%  Cl)  number  of  Trichoptera  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  67.  Mean  (95%  Cl)  number  of  sliredder  taxa  for  Soda  Butte  Creek,  1995. 


129 


10  %96  PUB  UB0LU 


o 

CO 


CD 


o 

"O 


-B 


E 

00 

E 

CO 

0 


E 


CO 


I 

(A 

bO 

_c 


eg 

E 


4- 


CO 


ID  CN  ID  T-  to 

• • • 

CM  'c-  O 

Bxe\  J0pp0jqs  p # 


E 


o 

u. 

To 

o. 

3 


• . V"'  • 


.><•'• . . yK- 


<v»  /"’  V T ^''  - '> 


;tfr  •■ ' 


"'  't;.r 


rt'. 


fVfl 


.:2^- 


► •', 


* i^  iiiL 


-*i- 


Figure  68.  Mean  (95%  Cl)  number  of  scraper  taxa  for  Soda  Butle  Creek,  1995. 
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Figure  69.  Mean  (95%  Cl)  number  of  collector-gatherer  and  filter-gatlierer  iaxa  for  Soda  Butte  Creek,  1995. 


131 


10  %96  PUB  UB0LU 


CO 


O^OOh-COlO'^COCNJ 

Bxe\  6]  puB  6o  jO  # 


upstream  Qlfp  downstream 


Figure  70.  Mean  (95%  Cl)  number  of  predator  taxa  for  Soda  Butte  Creek,  1995. 


132 


10  %96  PUB  UB9UJ 


BXBJ  JOlBpejd  p # 


downstream 


Figure  71.  Mean  (95%  Cl)  % shredder  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  72.  Mean  (95%  Cl)  % scraper  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  73.  Mean  (95%  Cl)  % collector-gatherer  and  filter-gatherer  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  74.  Mesn  (95%  Cl)  % predator  taxa  for  Soda  Dutte  Creek,  1995. 
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Figure  75.  Mean  (95%  Cl)  % dominant  taxa  for  Soda  Butte  Creek,  1995. 
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Figure  76.  Mean  (95%  Cl)  % EPT  taxa  for  Soda  Butte  Creek,  1995. 
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Table  33.  Mean  C±  95%  Cl),  variance,  standard  deviation  (s.d.),  and  coefficient  of 
variation  (CV)  for  metrics  calculated  for  seven  reference  streams  (N  = 34  samples, 
except  Metals  Tolerance  Index,  N = 17  samples) 


metric 

mean 

(±  95%  C.I.) 

variance 

s.  d. 

CV 

a of  organisms 

393.00 
C±  60.36) 

32232.73 

179.53 

0.46 

taxa  richness' 

25.59 

(±1.31) 

15.28 

3.91 

0.15 

EPT  richness' 

20.53 

(±0.91) 

7.35 

2.71 

0.13 

# of  Plecoptera 
taxa 

3.59 

(±0.45) 

1.76 

1.33 

0.37 

it 

Ephemeroptera 

taxa 

9.79 

(±  0.40) 

1.38 

1.17 

0.12 

it  of  Trichoptera 
taxa 

7.15 

(0.54) 

2.55 

1.60 

0.22 

a of  shredder  taxa 

1.85 

(±0.38) 

1.28 

1.13 

0.61 

it  of  scraper  taxa 

3.65 

(±0.45) 

1.81 

1.35 

0.37 

it  of  eg  and  fg 
taxa^ 

8.29 

(0.57) 

2.88 

1.70 

0.20 

it  of  predator  taxa 

6.18 

(±0.59) 

3.06 

1.75 

0.28 

% shredder  taxa 

6.86 

(±1.23) 

13.38 

3.66 

0.53 

% scraper  taxa 

14.27 
(±  1.58) 

22.02 

4.69 

0.33 

% eg  + fg  taxa^ 

33.10 

(±2.68) 

63.76 

7.98 

0.24 

% predator  taxa 

24.22 

(±1.99) 

35.33 

5.94 

0.24 

% dominant  taxa 

ifflP 

1.38 

(±0.14) 

0.18 

0.42 

0.30 

Shannon  Diversity 
Index 

2.51 

(±0.17) 

0.17 

0.41 

0.16 

%EPT  richness' 

80.81 
(±  2.54) 

0.57 

7.56 

0.09 

%EPT 

organisms 

86.98 

(±3.23) 

0.92 

9.62 

0.11 

% Chironomidae 

5.33 

(±  1.49) 

0.19 

4.43 

0.83 

Metals  Tolerance 
Index 

0.69 
(±  0.06) 

0.03 

0.18 

0.27 

'Taxa  richness  is  the  number  of  taxa  per  sample. 

^EPT  richness  is  the  number  of  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa  in  a sample, 
^cg  and  fg  is  the  number  of  collector-gatherer  and  filter-gatherer  taxa 
* HBI  is  Hilsenhoff  Biotic  Index 
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Table  34.  Mean,  variance,  standard  deviation  (s.d.),  coefficient  of  variation  (CV)  and 
95%  confidence  interval  (Cl)  for  metrics  calculated  for  Hyalite  Creek,  1994. 


mean 

variance 

s.  d. 

CV 

Cl 

Ephemeroptera  taxa 

8.00 

0.67 

0.82 

0.10 

0.80 

Plecoptera  taxa 

2.00 

0.67 

0.82 

0.41 

0.80 

Trichoptera  taxa 

7.50 

0.33 

0.58 

0.08 

0.57 

taxa  richness 

21.00 

11.33 

3.37 

0.16 

3.30 

organisms 

187.75 

667.58 

25.84 

0.14 

25.32 

shredder  taxa 

1.00 

1.33 

1.15 

1.15 

1.13 

scraper  taxa 

3.00 

0.67 

0.82 

0.27 

0.80 

eg  and  fg  taxa 

7.00 

0.67 

0.82 

0.12 

0.80 

predator  taxa 

3.50 

0.33 

0.58 

0.16 

0.57 

% dominant  taxa 

0.34 

0.01 

0.09 

0.25 

0.08 

% shredder  taxa 

0.04 

0.00 

0.05 

1.16 

0.05 

% scraper  taxa 

0.15 

0.00 

0.06 

0.43 

0.06 

%cg  and  fg  taxa 

0.34 

0.00 

0.06 

0.18 

0.06 

% predator  taxa 

0.17 

0.00 

0.05 

0.29 

0.05 

shannon  diversity  index 

2.20 

0.05 

0.23 

0.69 

0.23 

%EPT  taxa 

0.84 

0.01 

0.07 

0.09 

0.07 

EPT  taxa 

17.50 

3.00 

1.73 

0.10 

1.70 

% EPT  organisms 

0.97 

0.00 

0.01 

0.01 

0.01 

% Chirinomidae 

0.00 

0.00 

0.00 
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0.00 
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Table  35.  Mean,  variance,  standard  deviation  (s.d.),  coefficient  of  variation  (CV)  and 
95%  confidence  inter\^al  (Cl)  for  metrics  calculated  for  Squaw  Creek,  1995. 


mean 

variance 

s.d. 

CV 

Cl 

Taxa  richness 

23.00 

3.00 

1.73 

0.08 

1.96 

# organisms 

243.67 

10326.33 

101.62 

0.42 

115.00 

EPT  taxa 

20.67 

5.33 

2.31 

0.11 

2.61 

Ephemeroptera  taxa 

9.00 

1.00 

1.00 

0.11 

1.13 

Plecoptera  taxa 

3.67 

0.33 

0.58 

0.16 

0.66 

Trichoptera  taxa 

9.00 

1.00 

1.00 

0.11 

1.13 

“/oEphemeroptera  taxa 

39.00 

5.33 

2.31 

0.06 

2.61 

“/oPlecoptera  taxa 

16.00 

3.00 

1.73 

0.11 

1.96 

%Trichoptera  taxa 

34.67 

8.33 

1 

2.89 

0.08 

3.27 

%EPT 

89.67 

16.33 

4.04 

0.05 

4.57 

%dominant  taxa 

31.33 

25.30 

5.03 

0.16 

5.69 

“/oEPTindividuals 

96.50 

0.09 

0.31 

0.00 

0.35 

# shredder  taxa 

1.00 

0.00 

0.00 

0.00 

0.00 

% shredder  taxa 

0.04 

0.00 

0.00 

0.08 

0.00 

# scraper, shredder  taxa 

6.00 

0.00 

0.00 

0.00 

0.00 

% scraper, shredder  taxa 

26.19 

4.25 

2.05 

0.08 

2.32 

# eg  and  fg  taxa 

5.67 

2.33 

1.53 

0.27 

1.73 

% eg  and  fg  taxa 

24.41 

25.87 

5.09 

0.21 

5.76 

# scraper  taxa 

6.00 

0.00 

0.00 

0.00 

0.00 

% scraper  taxa 

26.19 

4.25 

2.06 

0.08 

2.33 

# predator  taxa 

7.33 

2.33 

1.53 

0.21 

1.73 

% predator  taxa 

31.75 

24.92 

4.99 

0.16 

5.65 

shannon  diversity  index 

2.37 

0.16 

0.16 

0.07 

0.18 

% Chironomidae 

0.02 

0.02 

0.02 

0.81 

0.02 
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Table  36.  Mean,  variance,  standard  deviation  (s.d.),  coefficient  of  variation  (CV)  and 
95%  confidence  interval  (Cl)  for  metrics  calculated  for  Cottonwood  Creek,  1995. 


mean 

variance 

s.d. 

CV. 

Cl 

Taxa  richness 

24.00 

27.00 

3.00 

0.11 

3.39 

#organisms 

533.00 

762.00 

124.23 

0.18 

140.59 

EPT  taxa 

23.00 

23.00 

1.15 

0.05 

1.31 

Ephemeroptera  taxa 

10.00 

11.00 

0.58 

0.06 

0.65 

Plecoptera  taxa 

4.00 

3.00 

1.00 

0.25 

1.13 

Trichoptera  taxa 

9.00 

9.00 

0.58 

0.06 

0.65 

%Ephemeroptera  taxa 

0.52 

0.62 

0.05 

0.09 

1.06 

%Plecoptera  taxa 

0.04 

0.02 

0.01 

0.24 

0.01 

%Trichoptera  taxa 

0.43 

0.24 

0.10 

0.32 

0.12 

%EPT  taxa 

0.96 

0.85 

0.07 

0.08 

0.08 

%dominant  taxa 

0.36 

0.24 

0.06 

0.21 

0.07 

%EPT  individuals 

0.92 

0.00 

0.06 

0.07 

0.07 

# shredder  taxa 

3.00 

0.00 

0.00 

0.00 

0.00 

% shredder  taxa 

0.07 

0.00 

0.01 

0.11 

0.01 

# scraper  and  shredder 
taxa 

4.00 

0.00 

0.00 

0.00 

0.00 

% scraper  and  shredder 
taxa 

0.15 

0.00 

0.02 

0.11 

0.02 

# eg  and  fg  taxa 

9.33 

0.33 

0.58 

0.06 

0.65 

%cg  and  fg  taxa 

0.34 

0.00 

0.04 

0.11 

0.04 

# scraper  taxa 

2.00 

0.00 

0.00 

0.00 

0.00 

% scraper  taxa 

0.07 

0.00 

0.01 

0.11 

0.01 

# predator  taxa 

9.00 

1.00 

1.00 

0.11 

1.13 

% predator  taxa 

0.33 

0.00 

0.04 

0.12 

0.04 

% chirinomidae 

0.01 

0.01 

0.01 

0.93 

0.01 

shannon  diversity  index 

2.28 

0.02 

0.02 

0.01 

0.02 
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DISCUSSION 

Multivariate  Analysis  vs.  a Multi-metric  Approach 

Principal  components  analysis  (PCA)  for  samples  collected  from  Bozeman 
Creek  (seasonal  variability),  the  three  Bear  Creeks  (within-site  and  between  stream 
variability),  and  Pebble  Creek  and  Soda  Butte  Creek  (within-stream  variability)  were 
based  on  individual  species  relative  abundances  grouped  into  sampling  units.  PCA 
suggested  seasonal  species  differences  in  Bozeman  Creek  and  between-stream  species 
abundance  similarities  in  the  three  Bear  Creeks.  PCA  also  suggested  differences  and 
similarity  of  species  abundances  between  2nd  order  and  3rd  order  reaches  in  Pebble 
Creek  and  between  Rosgen  stream  types. 

The  multi-metric  approach  on  the  other  hand,  did  not  rely  on  individual  species 
relative  abundances.  The  multi-metric  approach  emphasized  community  structure  and 
function.  Multimetric  emphasis  on  community  structure  and  function  is  based  on  stream 
ecosystem  theory  that  hypothesizes  that  individual  species  abundances  (and  their 
presence  or  absence)  naturally  vary  between  and  within  streams,  while  the  structure 
and  function  of  a streams'  community,  under  natural  conditions,  should  be  predictable 
between  stream  reaches  within  an  ecoregion  (Vannote  et  al.  1980,  Cummins  1974). 

PCA  analysis,  as  was  used  in  this  study,  was  therefore,  a more  useful  tool  for 
measuring  changes  in  species  composition  between  seasons  and  longitudinally,  while 
the  multi-metric  approach  was  more  useful  for  measuring  changes  in  the  structure  and 


function  of  stream  communities. 
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Within-site  Variability 

The  % Chironomidae  was  the  only  'least  impaired'  reference  site  metric  that 
had  a mean  CV  (0.43  +.  1.01)  large  enough  to  be  considered  'too  variable'  for  use  in 
mountain  stream  rapid  assessment.  All  of  the  other  metrics  had  acceptable  within-site 
CV  variability  (<  0.23).  Many  state  and  federal  water  quality  agencies  expressed 
concern  about  the  high  variability  of  the  number  of  organisms  metric  and  its  usefulness 
in  rapid  assessment  (Bukantis  1995,  Rosenberg  and  Resh  1993).  The  mean  CV  for 
reference  streams  in  this  study  for  the  number  of  organisms  metric  was  0.23  which  is 
'not  unusual'  for  field  studies  (Clarke  1994).  Therefore,  from  this  study,  it  appeared 
that  within-site  variability  of  the  number  of  organisms  for  mountain  reference  streams 
in  Montana  was  low  enough  to  allow  it  to  be  useful. 

Within-stream  Variability 

Some  metric  trends,  for  example  % Chironomidae  organisms,  and  PC  A results 
from  2nd  order  and  3rd  order  streams  sites  in  Pebble  Creek  (1995)  suggest  that  stream 
order  particularly,  crenal  (cold,  headwater  1st  to  2nd  order)  sites  differ  enough  in 
macroinvertebrate  communities  to  possibly  influence  water  quality  assessment.  Which 
is  not  surprising.  Headwater  reaches  in  mountain  streams  typically  have  different 
macroinvertebrate  communities  than  downstream  reaches  (Gustafson  1990,  Ward  1995; 
Ward  and  Kondratieff  1992). 

Ward  (1995)  reported  that  downstream  faunal  changes  in  alpine  streams  were 
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dependent  on  maximum  water  temperatures  with  a higher  relative  abundance  of 
Chironomidae  found  in  the  headwaters.  Temperature  data  collected  in  this  study  was  a 
single  reading  for  each  site  and,  therefore,  not  sufficient  to  attempt  temperature 
influence  analysis  on  the  macroinvertebrate  communities. 

Although  the  % Chironomidae  metric  showed  the  most  within-site  variability,  it 
could  be  a useful  metric  for  detecting  where  the  change  in  macroinvertebrate 
communities  from  crenal  (cold,  headwater)  communities  to  rhithral  (mid-reach) 
communities  occurs  in  unimpaired  streams.  For  example,  in  Pebble  Creek,  mean  % 
Chironomidae  at  Km  2 was  54.18%  (+.  14.08%)  then  decreased  at  Km  4 to  23.50%. 
While  further  downstream  at  3rd  order  sites  Km  6.5,  Km  9,  Km  13,  and  Km  16  the  % 
Chironomidae  ranged  from  only  4%  (±,2%)  to  11%C±.2%).  The  large  abundance  of 
Chironomidae  organisms  also  affected  other  metrics  such  as  the  lowered  %EPT  taxa. 

Gustafson  (1990)  reported  a decelerating  change  in  community  structure  in 
streams  in  the  Gallatin  River  drainage  which  he  correlated  to  a decelerating  change  in 
stream  order.  Gustafson  (1990)  reported  that  macroinvertebrate  changes  occur  more 
rapidly  along  small  streams  than  along  similar  distances  of  larger  streams. 

Pebble  Creek  2nd  order  stream  reach  occurred  from  1 kilometer  to 
approximately  5 kilometers  from  origin,  which  included  site  2 Km  (2670  masl)  and  site 
4 Km  (2600  masl)  for  a total  length  of  about  4 km.  The  3rd  order  stream  reach 
occurred  in  Pebble  Creek  from  approximately  5 kilometers  to  18  kilometers  with  its 
confluence  with  Soda  Butte  Creek  (sites  Km  6.5,  Km  9,  Km  13,  and  Km  16)  for  a total 
length  of  about  13  kilometers.  Therefore,  the  decelerating  rate  of  change  in  stream 
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order  and  its  associated  change  in  macroinvertebrate  communities  should  be  considered 
in  rapid  assessment  in  mountain  streams  of  Montana,  with  particular  reference  to 
higher  elevation  headwater  (crenal)  sections.  Gustafson  (1990)  recommended  stream 
order  as  the  best  measure  of  stream  size  because  it  is  the  most  widely  used  and  it 
describes  faunal  response  to  stream  size  better  than  other  measures  (ex.  link  magnitude, 
drainage  area). 

The  differences  in  stream  types  (B3  vs  C4)  in  Pebble  Creek  may  have  affected 
some  of  the  metrics,  for  example,the  number  of  organisms  metric  and  number  of  taxa 
(taxa  richness).  The  lower  number  of  organisms  collected  from  Km  13  and  Km  16  (B3 
type),  as  compared  with  Km  2,  Km  4,  Km  6.5  and  Km  9 (C4  type),  may  have  been 
due,  in  part,  to  larger  substrate  size  and  less  availability  of  interstitial  spaces  available 
to  macroinvertebrates  in  the  B3  sites.  The  C4  type  sites  were  comprised  of  loose 
cobbles  and  gravels  and  therefore,  provided  much  more  habitat  for  macroinvertebrates, 
thus  their  higher  abundances.  Likewise  for  taxa  richness.  Stream  type  with  an 
emphasis  on  substrate  size  and  embeddedness  should  also  be  incorporated  into  future 
assessment  programs. 


Seasonal  Variability 

The  increase  in  the  number  of  organisms  between  early  summer  1994  and  late 
summer  1995  in  Bozeman  Creek  is  most  likely  explained  by  the  growth  and 
development  of  macroinvertebrates  eggs  and  larvae  that  were  deposited  in  the  spring 
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and  early  summer.  These  macroinvertebrate  larvae  thus  became  abundant  enough  and 
large  enough  to  be  collected  in  the  1 mm  mesh  kick  net  through  the  summer  season. 

The  slight  decrease  in  the  number  of  organisms  from  late  summer  to  fall  was  due  to  the 
adult  emergence  of  rapid  summer  developing  larvae  and  their  subsequent  absence  in  the 
benthos.  The  mean  number  of  taxa  for  these  three  sample  times  also  reflects  the 
presence  of  many  rapidly  maturing  taxa  throughout  the  summer,  their  emergence  as 
adults  in  late  summer,  and  their  absence  in  fall  samples.  %EPT  taxa  also  reflects  this 
phenomenon. 

The  increase  in  the  metric  values  of  % shredders  from  early  summer  to  autumn 
is  a result  of  the  growth  and  development  of  shredder  taxa  that  are  typically  most 
abundant  in  the  autunm  through  winter  seasons.  Shredder  taxa  have  evolved  to  exploit 
the  input  of  organic  matter  in  the  form  of  increased  leaf  litter  in  the  autumn  from 
deciduous  riparian  flora. 

Predator  taxa  are  present  in  the  benthos  at  a consistent  relative  abundance 
throughout  the  year.  The  slight  increase  in  predators  in  winter  and  spring  can  be 
attributed  to  the  increase  in  Plecoptera  taxa,  which  along  with  Rhyacophilidae 
(Trichoptera)  taxa,  are  the  primary  macroinvertebrate  predators  found  in  mountain 
streams  in  Montana. 

Many  of  these  observed  seasonal  changes  in  metrics  can  be  easily  explained  by 
natural  changes  in  stream  macroinvertebrate  communities  as  predicted  by  the  River 
Continuum  Concept  (Vannote  et  al.  1980)  and  general  stream  ecosystem  theory 
(Cummins  1974).  Although  the  results  presented  here  are  not  startling,  the  importance 
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of  this  seasonal  study  was  that  the  methods  used  in  rapid  bioassessment  for  mountain 
streams  in  Montana  were  sensitive  to  these  changes.  Variability  in  metric  values  due  to 
natural  seasonal  changes  in  macroinvertebrate  communities  can  easily  be  avoided  and 
usually  are  by  sampling  in  only  one  designated  season  (usually  summer)  . 

If,  though,  the  objective  of  a rapid  assessment  was  to  document  the  impact  of 
riparian  changes  from  logging  associated  activities,  then  it  might  be  worthwhile  to 
collect  macroinvertebrate  samples  in  the  fall  or  winter  seasons  when  the  number  of 
shredder  taxa  is  normally  higher.  An  upstream/downstream  or  before-and-after 
approach  would  be  appropriate. 

Spring  collections,  before  runoff,  should  yield  the  most  information  because 
there  usually  are  more  taxa  (and  EPT  taxa)  to  be  found  in  the  benthos  at  that  time  of 
year  (Gustafson  1990)  although,  this  smdy  did  not  entirely  support  this.  For  example. 
Spring  1995  Bozeman  Creek  taxa  richness  means  were  the  highest  at  32.67  (+.  3.97) 
while  Summer  1995  samples  had  a mean  of  32  (+.1.13)  taxa  . 

Impaired  vs.  Least  -impaired  (Heference")  Streams 

Because  the  % Chironomidae  metric  was  the  most  variable  (mean  CV  = 0.43 
+.  1.01),  had  25th  to  75th  percentile  overlap  between  reference  and  impaired  sites,  and 
only  comprised  about  5%  of  the  number  of  organisms  per  reference  site  sample,  the  % 
Chironomidae  organism  metric  should  not  be  used  at  this  time  for  use  in  assessing 
water  quality  impairment  in  mountain  streams.  Chironomidae  larvae  are  very  small, 
difficult  to  identify  at  the  species  level,  and  very  little  is  known  about  their  life  history 
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and  ecology  in  mountain  streams  of  Montana,  which  also  makes  their  use  limited 
(Gustafson  1990). 

Table  37.  summarizes  the  longitudinal  and  seasonal  variability  of  19  metrics.  It 
also  summarizes  the  direction  of  change  in  a metric  due  to  impairment  and  whether  the 
metric  appears  to  be  useful  for  impairment  detection  based  on  overlap  between  the  25th 
to  75th  percentiles  of  the  seven  reference  streams  and  the  three  impaired  streams. 

Squaw  vs.  Cottonwood  Creek 

Squaw  Creek  samples  did  not  show  impairment  as  well  as  Soda  Butte  Creek 
(site  2),  Hyalite  Creek,  and  Divide  Creek.  Although,  this  did  not  reflect  the  full 
impacts  of  logging  associated  activities  in  Squaw  Creek.  The  removal  of  most  of 
Squaw  Creek's  woody  debris,  at  least  in  the  lower  5 km,  appeared  to  have  caused 
Squaw  Creek  to  downcut  its  streambed.  Due  to  this  downcutting,  high  quality 
riffle/pool  sequences  were  hard  to  fmd  in  the  lower  5 km  of  Squaw  Creek.  High 
quality  riffle/pool  sequences  at  the  study  site  section  of  Cottonwood  Creek  on  the  other 
hand,  were  numerous  and  there  was  no  apparent  past  removal  of  woody  debris  from 
the  stream  channel.  Habitat  quality  differences  between  Squaw  Creek  and  Cottonwood 
Creek  was  visually  apparent.  Therefore,  macroinvertebrate  metric  values  if  used  alone 
may  not  reflect  all  levels  of  impairment  and  habit  assessment  measurements  (ex. 
riffle/pool  sequences)  should  remain  an  integral  part  of  rapid  assessment. 
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Soda  Butte  Creek 

Mine  leachate  from  McLaren  tailings  apparently  had  an  effect  on  metric  values. 
Although  the  numbers  of  taxa  of  shredders,  collector-gatherer  and  filter-gatherer,  and 
predators  were  reduced  due  to  impairment  at  site  2,  their  relative  abundances  were  not 
affected.  On  the  other  hand,  the  number  of  scraper  taxa  at  site  2 in  Soda  Butte  Creek 
was  reduced  to  0 taxa  (Figure  8).  This  was  probably  due  to  the  accumulation  of  heavy 
metal  precipitate  in  the  form  of  ferric  hydroxide  (Chadwick  1974)  on  the  substrate 
surfaces.  Scraper  taxa  rely  on  algae,  bacteria,  and  organic  detritus  on  the  substrate  as 
their  food  source.  Whether  the  elimination  of  scraper  taxa  was  due  to  heavy  metal 
precipitate  covering  scraper  food  supply  or  direct  toxicity  of  the  precipitate  to  scraper 
taxa  is  unknown.  Due  to  their  low  abundance,  it  is  unknown  if  the  macroinvertebrates 
collected  from  Soda  Butte  Creek  at  site  2,  directly  below  the  tailings,  were  a resident 
population  or  were  only  present  due  to  drift. 
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Table  37.  Comparison  of  metrics  from  the  seven  reference  and  three  impaired  streams:  including  direction  of  metric  values  due  to 
impairment,  overlap  in  25th  to  75th  percentile,  longitudinal  and  seasonal  trends,  and  potential  of  metric  for  continued  use  in 
assessment  of  water  quality  in  mountain  streams. 
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Ephemeroptera,  Plecoptera,  and  Trichoptera 
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Based  on  Pebble  Creek  longitudinal  data 
Based  on  Bozeman  Creek  seasonal  data 
collecter-gatherer  and  filter-gatlierer 
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comments 

potential  use  to 
detect  changes 
between  headwater 
and  mid-reach 
streams 

designed  for 
organic  enrichment 
impairment 

potential  of 
metric  for  use  in 
detection  of 
impairment 

not  useful 

useful 

not  useful 
for  types  of 
impairment  found 
in  this  study 

not  useful 

not  useful 

seasonal  trend* 

yes 

yes 

yes 

lowest  in  spring 

na 

longitudinal 

trend' 

rapid  downstream 
decrease  from  2nd 
to  3rd  order 
reaches 

yes 

rapid  decrease 
from  2nd  to  3rd 
order  reaches 

variable 

increase  rapidly 
from  2nd  to  3rd 
order  reaches 

overlap  in  25th  to 
75th  percentile 

yes 

no 

yes 

yes 

yes 

impariment 
direction  of 
metric  value 

none 

decrease 

none 

increase 

none 

Metric 

% Chironomidae 
organisms 

Shannon  Diversity 
Index 
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OQ 
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% dominant  taxa 

cd 
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‘ Based  on  Bebbie  L.reek  longitudinal  data 
^ Based  on  Bozeman  Creek  seasonal  data 
^ Ephemeroptera,  Plecoptera,  and  Trichoptera 
* Hillsenhoff  Biotic  Index 
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CONCLUSION 

Within-site,  within-stream,  and  seasonal,  natural  variation  exists  in  the 
distribution  and  abundance  of  macroinvertebrates  in  mountain  streams  of  Montana. 

Due  to  this  natural  variabililty,  present  rapid  assessment  programs  are  forced  to  use 
wide  confidence  limits  for  their  reference  stream  metrics.  This  forfeits  sensitivity  and 
ability  to  detect  impairment.  Given  this  natural  variability,  most  of  the  metrics  used  in 
this  study  were  affected  by  severe  water  quality  impairment.  If  the  following 
recommendations  based  on  this  study  and  future  studies  are  considered  then  confidence 
limits  (ranges)  can  be  reduced  and  a more  accurate  assessment  of  impairment  can  be 
made: 

1)  Macroinvertebrate  collections  should  be  made  within  a two  month  time 
frame,  preferably  in  autumn  or  late  winter/early  spring.  This  will  avoid  the 
large  variability  associated  with  summer  adult  insect  emergence  and  seasonal 
variability.  If  summer  collections  must  be  made,  then  July  and  August  would 
be  the  prefered  months  which  will  reduce  variability  due  to  seasonal  changes. 

2)  Develop  separate  metric  regimes  for  1st  to  2nd  order,  cold,  headwater  stream 
segments  and  for  2nd,  3rd,  and  4th  order  mid-reach  stream  segments. 

3)  Develop  separate  metric  regimes  (suites  and  values)  for  different  types  of 
impairment  eg.  mining  activities,  logging  associated  activities,  channelization, 
etc. 


4)  Incorporate  more  multivariate  analysis,  such  as  Principal  Components 
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Analysis  or  Canonical  Correspondence  Analysis,  into  data  analysis.  These  are 
the  prefered  data  analysis  methods  used  in  water  quality  assessment  in  Europe 
and  can  easily  be  incorporated  into  the  multi-metric  approach  used  in  the  U.S. 
Some  additional  problems  in  a rapid  assessment  approach  to  mountain  stream 
impairment  still  need  to  be  addressed: 

1)  Functional  Feeding  Group  (FFG)  stams  of  many  mountain  stream  taxa  has 
yet  to  be  described  and  many  taxa,  if  not  most,  are  not  obligate  FFG  feeders. 
Many  taxa  can  and  will  feed  on  whatever  food  source  is  most  available.  Also, 
many  taxa  change  to  different  FFG's  throughout  their  larval  development.  The 
validity  of  the  use  of  FFG's  describe  mountain  stream  community  response  to 
impairment  should  be  further  considered. 

2)  This  study  only  compared  'severely'  impaired  streams  with  reference 
streams.  The  ability  of  metrics  and  rapid  assessment  to  detect  impairment  on 
'slightly'  impaired  streams  or  streams  of  varying  levels  of  impairment  needs  to 
be  further  smdied. 

3)  The  redundancy  of  some  of  the  metrics  used  in  this  study  was  not  addressed. 
PC  A would  be  a good  tool  for  future  metric  redundancy  analysis. 

4)  Some  analysis  and  metric  development  based  on  individual  taxa  as  "indicator 
species"  could  also  add  to  rapid  assessment  sensitivity  to  impairment. 

5)  Mountain  stream  health  should  be  a priority  because  mountain  streams  are 
the  headwaters  of  many  streams  and  rivers  in  Montana.  Whatever  impairment 
occurs  in  mountain  streams  affects  water  quality  downstream.  It  is  easier  to 
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protect  water  quality  before  impairment  occurs  than  to  try  to  remediate 
impairment  after  it  occurs. 

Finally,  future  research  efforts  should  focus  on  describing  the  namral 
similarities  and  differences  of  macroinvertebrate  communities  in  mountain  streams 
throughout  Montana  in  relation  to  the  physical  and  chemical  properties  of  streams  and 
their  drainages  on  a sub-ecoregional  or  watershed  level.  This  would  provide  for  a much 
needed  larger  reference  stream  data  base  which  then  can  be  used  for  comparison  with 


impaired  streams. 
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Table  38.  Species  list  for  Bozeman  Creek  seasonal  samples  (S  1,S2, S3, S4  = summer  1994  samples;  Al,A2,A3=autumn  1994 
samples;  W1,W2,W3  = winter  1995  samples;  Spl,Sp2,Sp3=  spring  1995  samples;  Sl*,S2*,S3*=summer  1995  samples). 
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Table  39.  Species  list  for  Soda  Butte  Creek  samples,  1995  (3  samples  per  site). 
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Table  40.  Species  list  for  Pebble  Creek  samples,  1995  (from  headwaters  Km  2 downstream  to  Km  16) 
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X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Tipulidae  1 

Tipulidae  2 
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Table  41.  Species  list  for  N.F.  Bear  Creek,  M.F.  Bear  Creek,  and  Bear  Creek  1994 
and  1995  (3  samples  from  each  stream  in  1995,  1 sample  from  each  stream  in  1994) 


N.F. Bear  95 

94 

M.F.  Bear  95 

94 

Bear  95 

94 

sample  if 

1 

2 

3 

1 

1 

2 

3 

1 

1 

2 

3 

1 

EPHEMEROPTERA 

Baetidae 

Baetis  bicaudatus 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

B.  tricaudatus 

Ephemerellidae 

Caudatella  hystrix 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Drunella  doddsi 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

D.  coloradensis 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

D.  grandis 

D.  spinifera 

Serratella  tibialis 

X 

X 

X 

X 

X 

X 

X 

Heptageniidae 

X 

Cinygmula  sp. 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Epeorus  albertae 

E.  deceptivus 

X 

X 

X 

X 

X 

X 

X 

X 

X 

E.  grandis 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

E.  longimanm 

X 

X 

X 

X 

Rhithrogena  sp. 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Siphlonuridae 

Ameletus  sp. 

X 

X 

X 

X 

X 

X 

PLECOPTERA 

Capniidae 

X 

Chloroperlidae 

X 

Suwallia  sp. 

X 

X 

X 

X 

Sweltsa  sp. 

X 
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Nemouridae 

Podmosta  sp. 

Zapada  cincitipes 

Z.  oregonensis 

Perlidae 

Doroneuria  theodora 

Hesperoperla  pacifica 

Perlodidae 

Isoperla  sp. 

I.  fulva 

X 

X 

I.  petersoni 

Kogotus  modesta 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Megarcys  sp. 

X 

X 

X 

X 

X 

X 

X 

X 

Pictitiella  expanse 

Skwala  sp. 

Peltoperlidae 

Yoroperla  brevis 

X 

Taeniopterygidae 

Taenionema  pacificum 

TRICHOPTERA 

Arctopsychidae 

Arctopsyche  grandis 

X 

X 

Parapsyche  elsis 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Brachyceniridae 

Brachycentnis  americanus 

B.  occidentalis 

Glossosomatidae 

Glossosoma  sp. 

X 

X 

X 

Philopotamidae 
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Polycentrepodidae 

Rhyacophilidae 

Rhyacophila  albertae 

R.  angelita  group 

X 

X 

X 

X 

X 

X 

X 

R.  coloradensis  group 

R.  hyalinata 

X 

X 

X 

X 

X 

X 

X 

X 

R.  narvae 

X 

X 

X 

X 

R.  oreta 

R.  vaccua 

X 

X 

X 

X 

R.  valuma  and  pellisa 

R.  vemna  and  brunnea  group 

X 

X 

X 

X 

X 

X 

R.  verrula 

X 

X 

X 

R.  vofixa 

X 

X 

Uenoidae 

Neothrema  sp. 

X 

X 

X 

X 

Dicosmoecus  sp. 

X 

X 

X 

X 

Oligophlebodes  sp. 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

unidentified  Trichoptera  1 

X 

X 

X 

unidentified  Trichoptera  2 

X 

X 

X 

COLEOPTERA 

Elmidae 

Clepelmis  sp. 

X 

Lara  sp. 

Narpus  sp. 

Adult 

DIPTERA 

Athericidae 

X 

X 

Blephariceridae 

X 

X 

Ceratopogonidae 

X 

X 

X 
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Culicoides  sp. 

Chironomidae 

Tanytarsus  sp. 

X 

X 

X 

X 

X 

X 

X 

Diamesa  sp. 

Pagastia  sp. 

Cricotopus  nostococladius 

X 

Eukiefferiella  sp. 

Tvetenia  sp. 

Deuterophlebiidae 

Empididae 

X 

X 

Psychodidae 

Simuliidae 

X 

X 

X 

X 

X 

X 

X 

Tabanidae 

Tipulidae  1 

X 

X 

Tipulidae  2 

HEMIPTERA 

X 

TURBELLARIA 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

OLIGOCHEATA 

HIRUDINEA 
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Table  42.  Species  list  for  Divide  Creek  1995  and  1994,  Cottonwood  Creek  1995  and 
Squaw  Creek  1995. 


Divide  1995 

94 

Cottonwood 

1995 

Squaw  1995 

sample  # 

1 

2 

3 

1 

1 

2 

3 

1 

2 

3 

EPHEMEROPTERA 

Baetidae 

Baetis  bicaudatus 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

B.  tricaudatus 

Ephemerellidae 

Caudatella  hystrix 

X 

X 

X 

X 

X 

X 

Drunella  doddsi 

X 

X 

X 

X 

X 

X 

X 

X 

X 

D.  coloradensis 

X 

X 

X 

X 

X 

X 

X 

X 

X 

D.  grandis 

X 

X 

X 

D.  spinifera 

X 

. 

Serratella  tibialis 

X 

X 

X 

X 

X 

X 

Heptageniidae 

Cinygmula  sp. 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Epeorus  albertae 

X 

X 

E.  deceptivus 

X 

X 

X 

X 

X 

X 

E.  grandis 

X 

X 

X 

X 

X 

X 

X 

E.  longimanus 

X 

X 

X 

X 

X 

Rhithrogena  sp. 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Siphlonuridae 

Ameletus  sp. 

PLECOPTERA 

Capniidae 

Chloroperlidae 

X 

X 

X 

X 
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Suwallia  sp. 

Sweltsa  sp. 

X 

X 

X 

X 

Nemouridae 

Podmosta  sp. 

Zapada  cincitipes 

X 

X 

X 

Z.  oregonensis 

Perlidae 

Doroneuria  theodora 

X 

X 

X 

X 

X 

X 

Hesperoperla  pacifica 

Perlodidae 

X 

Isoperla  sp. 

I.  fulva 

X 

X 

X 

X 

X 

X 

I.  petersoni 

X 

X 

Kogotus  modesta 

X 

Megarcys  sp. 

X 

X 

X 

Picntiella  expansa 

Skwala  sp. 

Peltoperlidae 

Yoroperla  brevis 

Taeniopterygidae 

Taerdonema  pacificum 

TRICHOPTERA 

Arctopsychidae 

Arctopsyche  grandis 

X 

X 

X 

X 

Parapsyche  elsis 

Brachycentridae 

Brachycentrus  americanus 

X 

X 

B.  occidemalis 

Glossosomatidae 
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Glossosoma  sp. 

X 

X 

X 

Philopoiamidae 

Polycentrepodidae 

Rhyacophilidae 

Rhyacophila  albertae 

X 

X 

X 

X 

X 

R.  angelita  group 

X 

X 

R.  coloradensis  group 

X 

X 

X 

R.  hyalinata 

X 

X 

X 

R.  narvae 

X 

X 

X 

R.  oreta 

R.  vaccua 

R.  valuma  and  pellisa 

X 

R.  vemna  and  brunnea  group 

X 

X 

X 

X 

X 

X 

R.  verrula 

X 

X 

X 

R.  vofixa 

X 

X 

Uenoidae 

Neothrema  sp. 

X 

X 

X 

X 

X 

Dicosmoecus  sp. 

Oligophlebodes  sp. 

X 

X 

X 

unidentified  Trichopiera  1 

X 

X 

X 

X 

X 

unidentified  Trichoptera  2 

X 

X 

COLEOPTERA 

Elmidae 

Clepelmis  sp. 

Lara  sp. 

Narpus  sp. 

Adult 

DIPTERA 

Athericidae 
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Blephariceridae 

- 

Ceratopogonidae 

Culicoides  sp. 

Chironomidae 

X 

X 

X 

Tanytarsus  sp. 

X 

X 

X 

X 

X 

Diamesa  sp. 

Pagastia  sp. 

Cricotopus  nostococladius 

X 

X 

X 

Eukiejferiella  sp. 

Tvetenia  sp. 

Deuterophlebiidae 

Empididae 

Psychodidae 

Simuliidae 

X 

X 

X 

X 

X 

Tabanidae 

Tipulidae  1 

X 

X 

X 

Tipulidae  2 

X 

HEMIPTERA 

TURBELLARIA 

OLIGOCHEATA 

X 

HIRUDINEA 
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